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Abstract 


Swarms  of  flying  robotic  insects  could  revolutionize  hazardous  environment  explo¬ 
ration,  search  and  rescue  missions,  and  military  applications.  Reducing  size  to 
insect  scale  enables  entrance  into  extremely  narrow  spaces  with  inherent  stealth 
advantages.  For  mass  production,  these  vehicles  must  have  reliable  and  repeatable 
fabrication  processes  that  define  flapping  wing  mechanisms  with  microscale  fea¬ 
tures  and  produce  large  flapping  amplitudes  at  frequencies  in  the  range  of  many 
insects.  This  report  focuses  on  the  design  and  fabrication  processes  of  flapping 
wing  mechanisms  for  these  types  of  robots. 

First,  the  design,  fabrication,  modeling,  and  experimental  validation  of  the 
Penn  State  Nano  Air  Vehicle  (PSNAV),  a  NAV  scale  piezoelectrically  actuated 
clapping  wing  mechanism,  is  presented.  A  flexure  hinge  allows  passive  wing  ro¬ 
tation  for  the  clapping  wing  mechanism.  Analytical  models  of  wing  flapping  and 
rotation  are  derived  and  validated  using  experimental  wing  trajectory  results.  The 
PSNAV  prototype  is  experimentally  shown  to  provide  approximately  54  deg.  peak 
to  peak  wing  rotation,  14  deg.  peak  to  peak  flapping  angle,  and  0.21  mN  of  thrust 
at  9.5  Hz.  At  25.5  Hz ,  the  prototype  produces  a  maximum  of  1.34  mN  of  thrust. 
The  PSNAV  model  accurately  predicts  the  wing  resonances  in  the  experimental 
prototype.  Model-predicted  thrust  is  lower  than  the  experimentally  measured  val¬ 
ues,  however. 


Towards  a  compliant  mechanism,  the  next  stage  of  this  research  introduces  a 
simple  process  to  monolithically  fabricate  flying  robotic  insects  at  the  pico  air  ve¬ 
hicle  (PAV)  scale  from  SUEX  dry  him,  an  epoxy  based  negative  photoresist  similar 
to  SU-8.  The  developed  process  has  fewer  steps  compared  to  other  methods,  does 
not  use  precious  metals,  and  greatly  reduces  processing  time  and  cost.  It  simul¬ 
taneously  defines  the  PAV  airframe,  compliant  happing  mechanism,  and  artihcial 
insect  wing  using  photolithography.  Using  this  process,  we  designed  and  fabricated 
the  LionFly,  a  happing  wing  prototype  actuated  by  a  PZT-5H  bimorph  actuator. 
Several  LionFly  prototypes  were  fabricated  and  experimentally  tested.  Theoret¬ 
ical  and  experimental  results  have  excellent  agreement  validating  the  compliant 
mechanism  kinematics  and  aerodynamic  added  mass  and  damping.  High  voltage 
tests  show  a  peak  to  peak  happing  angle  of  55  deg.  at  150  V  amplitude  with  150  V 
DC  offset  at  51  Hz  resonance.  Consistent  performance  from  multiple  prototypes 
demonstrate  the  reliable  and  repeatable  nature  of  the  fabrication  process. 

Lastly,  this  research  presents  detailed  modeling  and  experimental  testing  of 
wing  rotation  and  lift  in  the  LionFly.  A  flexure  hinge  along  the  span  of  the  wing 
allows  the  wing  to  rotate  in  addition  to  happing.  A  linear  vibrational  model  is 
developed  and  augmented  with  nonlinear  aerodynamic  forces  using  the  blade  cl¬ 
ement  method.  This  model  is  validated  using  experimental  testing  with  a  laser 
vibrometer  and  accurately  predicts  small  amplitude  wing  dynamics  in  air  and  vac¬ 
uum.  Strobe  photography  and  high  definition  image  processing  is  used  to  measure 
high  amplitude  wing  trajectories.  At  higher  amplitudes,  the  model  can  sufficiently 
predict  wing  trajectory  amplitudes,  but  phase  measurement  and  simulation  have 
slight  error.  The  LionFly  produces  46  deg.  hap  and  44  deg.  rotation  peak  to  peak 
with  relative  phase  of  12  deg.,  and  maximum  lift  of  71  fiN  at  37  Hz.  By  reducing 
the  inertia  of  the  wing  and  tuning  the  rotational  hinge  stiffness,  a  redesigned  device 


is  simulated  to  produce  lift  to  weight  ratio  of  one. 
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Chapter  1 
Introduction 


Miniature  unmanned  air  vehicles  have  been  increasingly  popular  over  the  last 
decade  for  potential  reconnaissance  and  search  and  rescue  applications.  DARPA 
defines  Micro  Air  Vehicles  (MAVs)  as  having  wingspans  less  than  15  cm,  and 
the  next  generation  Nano  Air  Vehicles  (NAVs)  as  having  wingspans  less  than  7.5 
cm,  weight  less  than  10  grams,  with  payload  capability  of  2  grams[2].  Recently, 
Wood  has  defined  a  Pico  Air  Vehicle  (PAV),  with  a  maximum  takeoff  weight  of 
500  milligrams  and  maximum  dimension  of  five  centimeters  [3].  Reducing  size  to 
insect  scale  enables  entrance  into  extremely  narrow  spaces  with  inherent  stealth 
advantages. 

Typical  air  vehicle  configurations  consist  of  fixed  wings[4]  or  rotary  wings [5]. 
The  superior  maneuverability  of  small  birds  and  insects  have  motivated  researchers 
to  explore  flapping  wing  mechanisms  for  NAVs [6,  7,  8].  The  potential  ability  to 
fly  at  both  high  and  low  speeds,  hover  and  glide  give  flapping  wing  air  vehicles 
advantages  over  miniature  fixed  wing  aircraft.  In  addition,  if  designed  properly, 
these  vehicles  could  hide  in  plain  sight  and  mimic  a  small  bird  or  insect  which 
current  rotary  wing  aircraft  cannot  achieve.  The  smaller  dimensions  and  unique 
wing  trajectory  specifications  present  many  challenges  in  the  design  and  fabrication 
of  flapping  wing  actuators  and  mechanisms. 
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Flying  insects  have  been  a  crucial  source  of  inspiration  in  the  design  of  flapping 
mechanisms  [9,  8,  10,  7].  Dipteran  insects  use  indirect  flight  muscles  that  use  the 
thorax  to  amplify  flapping  wing  motion  and  operate  at  the  natural  frequency  of 
the  combined  system  [8] .  In  addition  to  flapping  angles,  the  wings  undergo  prona¬ 
tion  and  supination,  or  rotation,  which  is  key  for  lift  production  [11],  Swarms  of 
these  flying  robotic  insects  could  revolutionize  hazardous  environment  exploration, 
search  and  rescue  missions,  and  military  applications.  For  mass  production,  these 
vehicles  must  have  reliable  and  repeatable  fabrication  processes  that  define  flapping 
wing  mechanisms  with  microscale  features  and  produce  large  flapping  amplitudes 
(>  40°)  at  frequencies  in  the  range  of  many  insects  (35-100  Hz)  [12]. 

This  research  focuses  on  the  subsystem  of  flapping  wing  mechanisms,  including 
the  actuator,  amplification  mechanism,  and  wing;  the  interested  reader  is  directed 
to  reviews  of  state-of-the-art  energy  storage,  power  electronics,  sensing,  and  other 
payload  requirements  for  NAV/PAVs  [13,  14,  15].  The  following  section  presents 
background  and  comparison  on  the  state  of  the  art  in  flapping  wing  actuators  and 
mechanisms. 


1.1  Background  of  Flapping  Wing  Mechanisms 

1.1.1  Literature  Review 

At  the  MAV  scale,  many  fabrication  and  actuation  methods  are  similar  to  the  Delfly 
Micro  from  the  Delft  Univerisity  of  Technology.  This  10  cm  wingspan  MAV  is  hand 
built  using  carbon  rods,  balsa  wood,  Mylar  foil,  and  powered  by  a  DC  electric  motor 
coupled  to  a  gearbox  and  linkage  mechanism  using  conventional  pin  joints[16]. 
Coupled  to  a  MAV  scale  motor  and  gearbox,  Pornsin-Sirirak  et  ah  microfabricate 
7  cm  Ti-alloy  and  parylcne  C  wings  using  wet  etching[17].  AeroEnvironment  has 
recently  fabricated  the  “Nano  Hummingbird”  which  is  capable  of  sustained  flight 
and  reconnaissance  [18]. 
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Wood  et  al.  [19]  have  developed  the  Smart  Composite  Microstructure  (SCM) 
process  where  carbon  fiber  and  polyimide  layers  are  laser  micromachined  and  folded 
to  form  compliant  mechanisms,  enabling  fabrication  and  assembly  of  the  Harvard 
Microrobotic  Fly  (HMF),  a  60  mg,  3-crn  wingspan  PAV  that  uses  a  compliant 
flapping  wing  transmission  to  create  large  wing  angles  from  piezoelectric  bimorph 
displacement [8].  Successful  tethered  and  externally  powered  flight  of  the  HMF 
has  shown  that  piezoelectric  actuation  with  compliant  mechanism  amplification  is 
suitable  for  PAV  flight.  Using  MEMS  fabrication,  Bronson  et  al.  [20]  fabricate  the 
smallest  NAV  in  the  literature,  using  PZT  thin  film  actuators  directly  coupled  to  a 
2.5  mm  Si02/ Si3N4/Ti-Au  wing  that  produces  large  flapping  angle  at  resonance. 
Bao[21]  et  al.  microfabricated  a  electromagnetically  powered  flapping  wing  NAV 
from  multiple  layers  of  SU-8  photoresist. 

Many  actuation  methods  exist  in  the  design  space  of  flapping  wing  MAVs/NAVs: 
rubber  bands  [22],  reciprocal  chemical  muscles [23],  shape  memory  alloys [24,  25], 
electromagnetic  actuators [10],  electromagnetic  motors[6,  7,  26,  27,  28]  and  piezo¬ 
electric  actuators[8,  29,  30,  31].  Electromagnetic  motors  have  been  the  popular 
choice  for  larger  scale  mechanisms  [17,  9].  For  PAVs,  linear  electromagnetic  ac¬ 
tuation  [21]  and  bulk  PZT  bimorph  actuators  [8],  and  thin  film  PZT  unimorph 
actuators  [20]  have  been  used. 

Wing  rotation  is  crucial  for  insect  flight  and  has  been  extensively  studied  by 
biologists  and  engineers  [11,  32,  33,  34],  To  produce  aerodynamic  forces  sufficient 
for  flight,  insects  have  large  flapping  angles  (>  40°),  high  flapping  frequencies 
(>  40 Hz),  and  wing  twist,  or  rotation,  at  proper  times  during  the  flapping  cycle 
[12].  Many  insects  do  not  use  significant  muscle  power  in  rotating  their  wings,  but 
aerodynamic  and  inertial  forces  cause  passive  wing  rotation  [11].  Many  compu¬ 
tational  fluid  dynamic  studies  have  been  conducted  regarding  the  complex  fluid 
structure  interaction  on  wing  rotation  and  flexibility  [35,  36].  In  miniature  air 
vehicle  design,  wing  rotation  has  been  achieved  by  direct  actuation  [29],  elastody- 
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namic  tuning  [30],  and  aerodynamic  loading  [8].  From  a  vehicle  design  standpoint, 
achieving  wing  rotation  passively  is  preferred  for  reduced  system  mass,  complex¬ 
ity,  and  power  load.  Designing  this  type  of  vehicle  requires  accurate  modeling  of 
aerodynamic  and  inertial  loads  on  the  wing  during  the  flapping  cycle. 

Recently,  engineers  have  used  quasi-steady  aerodynamic  models  using  the  blade 
element  method  to  elucidate  aerodynamic  forces  and  moments  acting  upon  the 
wing  [33].  Nonlinear  analytical  models  have  been  created  which  also  describe 
flapping  and  rotation  dynamics  including  NAV  flapping  mechanisms  and  actuators 
[37].  Although  the  NAV  system  dynamics  are  nonlinear,  linear  systems  analysis 
allow  a  wide  variety  of  design  tools  and  have  shown  to  be  accurate  over  a  large 
operating  region  [38]. 

1.1.2  Component  Selection 

The  insect  flight  apparatus  consists  of  three  major  components:  the  flight  muscle, 
thorax,  and  wing.  Flight  muscles  in  insects  such  as  flies,  bees,  wasps,  and  beetles 
deform  the  thorax  amplifying  the  limited  muscle  contraction  into  large  flapping 
angles  [39].  For  pico  air  vehicle  design,  the  selection  of  the  actuator,  or  flight 
muscle,  determines  the  type  of  amplification  mechanism,  or  thorax,  required  for 
the  most  effective  flapping  wing  mechanism. 

Flapping  wing  mechanism  actuators  provide  either  linear  or  angular  displace¬ 
ment  which  must  be  converted  to  large  flapping  angles.  Electric  motors  are  most 
commonly  used  for  rotary  motion  input  with  gearbox  and  four  bar  mechanism 
amplificiation  using  conventional  pin  joints.  At  the  PAV  scale,  the  gearbox  be¬ 
comes  difficult  to  fabricate  in  an  integrated  way,  and  the  pin  joints  show  increased 
friction  losses.  Linear  actuators  have  been  used  with  conventional  pin  joints  [31], 
compliant  mechanisms  [30,  40,  8],  and  with  no  amplification  mechanism  [20,  41]. 
Devices  which  do  not  use  an  amplification  mechanism  have  actuators  which  attach 
to  the  wing  itself,  and  use  the  length  to  amplify  tip  displacement.  These  type  of 
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actuators  do  not  show  promise  because  no  similar  type  of  mechanism  is  found  in 
biology,  and  no  fabricated  device  has  been  shown  to  produce  lift  [20,  41].  On  the 
other  hand,  reduced  wear,  friction,  and  backlash,  compatibility  with  MEMS  pro¬ 
cessing  techniques  and  proven  PAV  application  [8]  motivate  the  use  of  compliant 
mechanisms. 


CM 


O 

cfa 

xs 

CD 

4^ 

Oh 

cd 

TJ 

cd 


fH 

CD 


f-{ 

O 

4^ 

cd 

2 

4^ 

cd 

cd 

0) 

3 


cd 

3 

xS 

Pi 

cd 

CD 


n 

o 

j} 

cd 

a 


O 

CD 


£ 


cd 

c? 


XI 

£ 


C/3 

CD 

4P> 

o 


CD 

£ 

.2 

O 

£6 

H 


-4H 

T3 


T3 

Pi 

cd 

m 


p 

.2 

o 

qj 

qp 

qj 

Q 


OJ 

O 

fH 

3 


CP 

a 

cd 

X 

H 


oj 

O, 

Eh 

?H 

o 

4— 1 

Cd 

Pi 

4^ 

CD 


§1 

Co  |  ~i 

43  r—H 

!> 


!> 


bO 


bO 


p?  cd 

g  fn 

l| 

O  ° 
£  ^ 
3  S 


K>~5 

4-1 

■a  & 

ch  ojO 
Pi  rrt 

oj  “ 
^  o 

in3  buO 

QJ  rH 

X  ‘  7^ 

®  P 

(~H  ?H 

Kfl  p> 

•SP  a 

-P  O 

K*^  fX 

£  .SP 

>  K 


CD 

b€ 

cd 

Jl 

o  bO 

QJ  <H 

'o  "43 

w  Id 

—I  fH 

cd  QJ 

■c  4 

gj  O 

P  PP 

P  £0 

OJ 

bO  >> 

fH  fH 

cd  .  QJ 

hP  > 


-P 

bO 


-P 

bO 


PP 

bO 


qj  jrj  41  K  hP  ffi 


-p  -P 

bO  bO 


-P 

bO 


ffi  ffi  4  ffi 


TP 

QJ 


a  a  m 

p  >  p  p  p  ■ p 

pp  P  bJO  bO  pd  — 

aj  i-P  ffi  K  p 

S  s  $ 


a  w)  a  ao  a 

.2  £  2p  .2  2p  .2 

QJ  hP  fD  OJ  fH J  QJ 
.OJ  5  _QJ  ferH 
^  > 


£ 


P=T  00 

fH  QJ 

2  qp 


-p  , _ , 

.CM,  ^ 
-p 

£  TT 

•  Dh  <D 


£ 

CD 

4— 1 

X 

CD 

> 


^H 

cd  4D 


o 

CD 


JD 

o 

C/3 

3 

a 

"cd 

*o 

eg 

4P> 

cd 


o 


K^5 

o 

< 

^~5 

f-H 

O 


CD 

£ 

b€  ^  a 

cd  cd  a; 

a^s 

o  o 

?H  ?-H 


4h 

CD 


CO 

cd 

H 


o  o 

CD  CD 


^  cd 
CD 

o 


CD 

Oh  X  *H 
Cd  S  CD 

W  W  co  P-,  H  Q 


CD 


6 


fH 

4H 

O 

r2 

'qj 

o 

N 

QJ 


p 

o 

tn 

'p 

cd 

& 

a 

o 

U 

bi 


^QJ 

3 


O1 


Pb 


o 

a. 


~e 


aj 

Pb 

a. 


p 


o 

o  ^ 

O  10  H 

p  a  cd 

r —  OJ  CD 

co  co  a 
odd 


CM  lO 
CM  QO 
H  CM 


10 


000 

O  O  H 
to  10  CO 

t-  00 


00  M4  LO 
QO  T — I 


-p 

to 

T— 1 

'  p . 

t-H 

ffi 

Eh 

-p 

Ph 

H 

H 

£ 

IS] 

M 

CS3 

CP 

Plh 

Plh 

PMN-PT  [46,  47]  14  8060  -1330  0.73  35-44 

PVDF  [46,  48]  9  1780  20  0.17  13-25 


7 


Many  types  of  linear  actuators  exist  that  are  compatible  with  compliant  ampli¬ 
fication  mechanisms.  To  reduce  the  scope  to  microrobotics,  six  candidate  actuator 
technologies  were  selected  for  comparison:  electromagnetic  [10],  electrostatic  [43], 
shape  memory  alloy  (SMA)  [24],  piezoelectric  [44],  thermal  [41],  and  dielectric  elas¬ 
tomers  [45].  Table  1.1.1  shows  a  qualitative  comparison  of  these  actuators  based  on 
the  review  and  analysis  by  Karpelson  and  Wood  [42] .  The  results  of  this  compari¬ 
son  may  change  based  on  the  specific  design  of  individual  actuators,  however  this 
is  useful  in  finding  fundamental  shortcomings  and  advantages  to  different  actua¬ 
tor  methodologies.  For  PAV  flapping  mechanisms,  electrostatic  actuators  do  not 
provide  the  requisite  force  and  deflection,  but  they  may  be  feasible  for  10  mg  scale 
microrobots  [49].  Electromagnetic  linear  actuators  are  strong  candidates,  however 
at  PAV  geometries  and  below  resistive  heat  losses  dominate  electromagnetic  force 
generation,  reducing  efficiency  [50].  SMA  actuators  have  very  high  energy  density 
and  thermal  actuators  are  very  powerful,  but  neither  can  provide  the  high  flapping 
frequencies  required  for  small  insect  flight.  Additionally,  the  mode  of  operation  of 
these  actuators  requires  Joule  heating  through  electrical  current  which  has  reduced 
efficiency  compared  with  other  methods.  Electromagnetic,  SMA,  and  thermal  ac¬ 
tuators  are  better  suited  for  larger  scale  applications  with  larger  dimensions  and 
lower  flapping  frequencies.  Dielectric  elastomers  are  very  promising  actuators  for 
PAV  applications,  however  require  voltages  an  order  of  magnitude  higher  than 
piezoelectric  actuators  which  already  suffer  from  high  voltage  requirements.  Re¬ 
cent  development  has  designed  and  fabricated  PAV  scale  power  electronics  suitable 
for  piezoelectric  actuators  [51].  The  piezoelectric  effect  enjoys  many  benefits  such 
as  favorable  scaling,  sufficient  force  and  deflection,  high  bandwidth,  and  good  effi¬ 
ciency  at  resonance  which  can  be  further  increased  with  charge  recovery  circuitry 

[52]. 

Table  1.1.2  shows  a  number  of  piezoelectric  actuator  materials  consisting  of 
piezoelectric  ceramics,  PZT-4  and  PZT-5H,  single  crystal  relaxor  ferroelectrics, 


PZN-PT  and  PMN-PT,  and  piezoelectric  polymer  PVDF.  Young’s  modulus,  E, 
density,  p,  piezoelectric  coefficient,  i,  electromechanical  coupling  coefficient, 
and  mechanical  quality  factor,  Qm  are  used  to  compare  these  materials.  Actua¬ 
tor  performance  is  dependent  on  configuration  such  as  unimorph,  bimorph,  rolled 
actuator  or  multilayer  stack,  but  can  be  fundamentally  compared  using  material 
properties.  Single  crystal  relaxor  ferroelectric  materials  have  impressive  piezoelec¬ 
tric  coefficient  and  electromechanical  coupling  coefficient,  however  poor  fracture 
properties  have  made  implementation  difficult  [44],  PVDF  has  very  small  piezo¬ 
electric  coefficient,  requiring  very  high  voltages  to  achieve  similar  strains  as  other 
piezoelectric  materials.  Hard  piezoelectric  ceramics  such  as  PZT-4  have  higher 
mechanical  quality  factor,  linear  response,  and  are  better  suited  for  resonance  ap¬ 
plications.  PZT-5H,  a  soft  piezoelectric  ceramic,  has  higher  piezoelectric  coefficient 
and  electromechanical  coupling,  but  large  hysteresis  and  nonlinearity.  Maximiz¬ 
ing  power  output  per  unit  weight  is  the  first  figure  of  merit,  making  PZT-5H  the 
best  piezoelectric  material  for  PAV  actuators.  For  an  additional  detailed  review  of 
piezoelectric  materials  and  actuators,  see  Appendix  A. 

Unimorph  [53],  bimorph  [44],  and  flextensional  [54]  actuators  have  been  used  for 
PAV  applications.  Unimorph  actuators  consist  of  an  active  piezoelectric  layer  (in 
the  case  of  a  bimorph  two  active  layers)  and  a  passive  elastic  layer  which  constrains 
the  field  induced  extension  of  the  active  layer  causing  bending  that  is  amplified  by 
the  length  of  the  beam,  increasing  tip  deflection.  Flextensional  devices  include  a 
compliant  mechanism  to  amplify  extension  of  the  piezoelectric  element  based  on 
the  geometry  of  the  mechanism.  None  of  these  actuators  provide  enough  deflection 
to  flap  a  wing  on  their  own,  and  they  must  have  additional  amplification  for  use 
in  PAV  mechanisms.  The  highest  force  and  deflection  output  per  unit  weight  is 
achieved  using  piezoelectric  bimorph  actuators,  especially  optimized  trapezoidal 
bimorphs  with  extensions  such  as  in  [44], 

The  linear  displacement  of  a  piezoelectric  bimorph  such  as  [44]  is  below  500  /im 
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amplitude  and  requires  significant  amplification  to  achieve  insect  like  flapping  am¬ 
plitudes  of  approximately  45  deg.  peak  to  peak  or  higher.  Flextensional  methods 
alone  provide  little  amplification  as  stated  before,  whereas  compliant  mechanisms 
based  on  traditional  four  bar  rigid  link  mechanisms  can  be  designed  to  provide 
sufficient  amplification  for  flapping  using  piezoelectric  actuation.  Wood  [55]  and 
Fearing  [40]  have  designed  three  dimensional  compliant  mechanisms  with  large 
amplification.  Designing  the  link  lengths  and  amplification  is  a  factor  of  the  indi¬ 
vidual  actuator  used  and  fabrication  constraints.  Jacobsen  [56]  discusses  in-plane 
mechanisms  which  create  out-of-plane  motion.  Creating  mechanisms  with  high 
amplification  in  this  way  can  be  significantly  more  difficult,  and  will  be  discussed 
in  this  report. 

1.2  Current  Research 

This  research  begins  with  the  design,  fabrication,  modeling,  and  experimental 
testing  of  the  first  NAV  scale  clapping  wing  mechanism  powered  by  piezoelectric  T- 
beam  actuators.  The  clapping  design  is  motivated  by  an  aerodynamic  mechanism 
called  the  Weis-Fogh  ’clap  and  fling’  utilized  by  certain  butterflies  and  wasps,  in 
which  opposing  wings  almost  touch  during  part  of  the  flap  cycle,  spawning  vortex 
structures  that  increase  thrust.  Many  different  clapping  wing  vehicle  designs  exist 
[26,  9,  7,  28,  57],  but  all  four- winged  ’X’  type  clapping  designs  use  a  mechanism 
consisting  of  electromagnetic  motors  with  gearboxes.  The  clapping  mechanism 
outlined  in  this  work  amplifies  the  T-bearn  displacement  using  miniature  revolute 
joints  and  hinges  spaced  a  certain  distance  apart,  creating  a  lever  mechanism.  The 
static  and  dynamic  operation  of  flapping  and  rotation  is  characterized  using  stobe 
photography  and  video,  and  lift  is  measured  using  a  force  transducer  coupled  with 
a  specially  designed  apparatus. 

Many  piezoelcctrically  actuated  flapping  wing  mechanisms  have  been  designed 
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and  fabricated  that  are  hand  assembled  and  have  3D  components  [8,  30,  31,  29]. 
Most  MEMS  compliant  mechanisms  operate  in-plane  with  in-plane  actuators  [58, 
59].  Jacobsen  et  al.  [56]  present  designs  for  lamina  emergent  mechanisms  that 
are  planar  mechanisms  capable  of  out  of  plane  motion.  This  research  introduces  a 
simple  flapping  wing  mechanism  fabricated  monolithically  from  SUEX  dry  him,  an 
epoxy  based  negative  photoresist  similar  to  SU-8.  The  developed  fabrication  pro¬ 
cess  has  fewer  steps  compared  to  other  methods,  does  not  use  precious  metals,  and 
greatly  reduces  processing  time  and  cost.  It  simultaneously  defines  the  PAV  air¬ 
frame,  compliant  happing  mechanism,  and  artihcial  insect  wing  using  photolithog¬ 
raphy.  Using  this  process,  we  designed  and  fabricated  the  LionFly,  a  happing 
wing  prototype  actuated  by  a  PZT-5H  bimorph  actuator.  Rapid  prototypes  were 
fabricated  with  precisely  dehned  features  and  material  properties  and  geometry 
that  are  similar  to  insects.  We  present  the  fabrication  process  and  characteriza¬ 
tion  of  SUEX,  static  and  dynamic  modeling  of  the  PAV,  and  experimental  results 
from  multiple  prototype  devices  that  demonstrate  applicability  of  the  developed 
fabrication  process  for  happing  wing  air  vehicle  applications. 

Finally,  in  miniature  air  vehicle  design,  wing  rotation  has  been  achieved  by 
direct  actuation  [29]  and  elastodynarnic  [30]  and  aerodynamic  [8]  passive  tuning. 
From  a  vehicle  design  standpoint,  passive  wing  rotation  is  preferred  for  reduced 
system  mass,  complexity,  and  power  consumption.  Designing  a  vehicle  with  passive 
wing  rotation  requires  accurate  modeling,  including  the  complex  interplay  between 
the  actuator,  mechanism  kinematics,  and  wing  vibration  and  aerodynamics  dur¬ 
ing  the  happing  cycle.  This  research  presents  detailed  modeling  and  experimental 
testing  of  wing  rotation  and  lift  in  the  LionFly.  The  goal  of  this  chapter  is  to  un¬ 
derstand  the  happing  and  rotation  dynamics,  and  the  lift-producing  mechanisms  in 
this  device.  A  linear  vibrational  model  is  developed  and  augmented  with  nonlinear 
aerodynamic  forces  using  the  blade  element  method.  Experimental  testing  using 
a  laser  vibrometer  in  air  and  in  vacuum  is  used  to  characterize  small  amplitude 
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flapping  and  rotation.  Strobe  photography  and  high  definition  image  processing 
is  used  to  measure  high  amplitude  wing  trajectories.  A  lift  measurement  system 
using  a  force  transducer  is  designed  and  used  to  measure  average  lift.  Finally,  the 
development  of  an  improved  fabrication  process  that  enables  thinner  flexures  and 
wings  is  motivated  by  simulations  using  the  experimentally  validated  model  that 
predict  a  lift  to  weight  ratio  of  one. 

1.3  Contribution  and  Organization  of  Report 

The  following  are  novel  contributions  in  the  field  of  piezoelectric  actuators,  flapping 
wing  mechanism  fabrication,  and  microfabrication: 

1.  Fabrication  of  piezoelectric  T-beam  actuators 

A  novel  fabrication  process  using  a  high  precision  dicing  saw  has  been 
developed  to  fabricate  piezoelectric  T-beam  acuators  at  the  millimeter  scale. 
The  static  and  dynamic  displacement  and  blocking  force  of  PZT-5H  piezo¬ 
electric  T-beam  acuators  have  been  modeled  and  experimentally  validated. 
The  understanding  developed  in  this  research  can  allow  piezoelectric  T-beam 
actuators  to  be  used  in  place  of  unimorph  or  bimorph  actuators  in  certain 
applications. 

2.  Clapping  wing  mechanism  for  NAV  applications 

Using  T-beam  actuators,  a  novel  clapping  wing  mechansim  has  been 
designed  and  fabricated  using  miniature  revolute  joints  made  from  diced 
medical  grade  syringes  and  tungsten  pins,  and  chemically  etched  stainless 
steel.  This  research  demonstrates  the  first  piezoelectrically  actuated  clapping 
wing  vehicle  to  produce  lift. 

3.  Monolithic  multilevel  fabrication  of  SUEX  dry  him 
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A  novel  fabrication  process  has  been  developed  to  monolithic-ally  fabri¬ 
cate  SUEX  dry  film.  The  process  exploits  the  optical  transmission  properties 
of  SUEX  dry  film,  defining  thin  flexible  regions  using  310  nm  wavelength  and 
thick  rigid  regions  using  365  nm  wavelength  UV  energy.  Low  cost  and  rapid 
prototypes  of  compliant  mechanisms  of  arbitrary  2D  dimensions  and  varied 
thickness  can  now  be  fabricated  using  photolithography. 

4.  The  LionFly,  a  SUEX  flapping  wing  mechanism 

A  novel  compliant  slider  crank  mechanism,  fabricated  monolithically 
from  SUEX,  has  been  developed  to  provide  large  flapping  and  rotation  an¬ 
gles.  Detailed  linear  vibrational  modeling  of  the  flapping  and  rotation  dy¬ 
namics  are  developed  and  are  augmented  by  nonlinear  aerodynamic  forces. 
The  model  is  experimental  validated  using  laser  vibrometer  measurement, 
strobe  photography  and  image  processing,  and  lift  measurement  with  a  force 
transducer  and  apparatus. 

The  report  is  organized  as  follows:  Chapter  2  describes  the  design,  fabrication, 
modeling  and  experimental  testing  of  the  clapping  wing  mechanism  powered  by 
piezoelectric  T-beam  actuators.  Chapter  3  describes  the  monolithic  fabrication  of 
SUEX  for  use  in  compliant  mechanisms,  and  the  design,  fabrication,  modeling,  and 
experimental  testing  of  a  compliant  slider  crank  mechanism  is  presented.  Chapter 
4  presents  detailed  modeling  and  experimental  testing  of  flapping  and  rotation 
in  the  LionFly  mechanism.  Chapter  5  concludes  with  results  from  the  different 
flapping  wing  mechanism  analyses. 


Chapter  2 

Clapping  Wing  Nano  Air  Vehicle 
Actuated  By  Piezoelectric 
T-b  earns 

2.1  Introduction 

In  this  chapter,  the  design,  fabrication,  modeling,  and  experimental  validation 
of  an  NAV  scale  clapping  wing  mechanism  is  presented.  Figure  2.1.1  shows  a 
conceptual  drawing  of  the  proposed  four  winged,  clapping  mechanism  for  an  NAV. 
A  flexure  hinge  allows  passive  wing  rotation  for  the  clapping  wing  mechanism.  To 
properly  predict  and  design  wing  trajectories,  analytic  models  of  wing  flapping  and 
rotation  are  derived,  including  unsteady  aerodynamic  wing  loading  corresponding 
to  the  low  Reynolds  numbers  flow.  The  simulated  NAV  wing  trajectories  and 
thrust  are  validated  using  the  PSNAV  prototype. 
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Figure  2.1.1.  Conceptual  drawing  (left)  of  the  PSNAV,  a  four  winged  piezoelectrically 
driven  nano  air  vehicle  and  fabricated  device  (right) 


2.2  Design  and  Fabrication 

2.2.1  T-Beam  Actuators 

Unimorph,  bimorph,  and  multilayer  actuators  have  been  developed  to  amplify 
PZT  strain  to  produce  larger  transverse  deflection  in  a  cantilevered  beam.  These 
actuators  require  process  steps  that  bond  separated  layers  of  piezoelectric  or  passive 
material  that,  in  bulk  piezoelectric  form,  are  difficult  to  translate  to  the  small  scale 
required  of  the  NAV  application.  T-beam  actuators  are  monolithically  fabricated 
from  a  bulk  piezoelectric  chip  and  provide  in-plane  and  out-of-plane  bending  in  a 
cantilevered  beam  [60].  T-Beam  actuators  have  T-shaped  cross  sections  and,  for 
this  application,  electrodes  on  the  top  and  bottom  of  the  structure.  Voltage  applied 
between  the  top  and  bottom  electrodes  produces  an  electric  field  that  concentrates 
in  the  web.  The  flange  material  is  essentially  passive  and  resists  contraction. 
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□  PZT 


Cr/Au 


(c)  ,-M— ,  , -  - !  |— M— | 

Figure  2.2.1.  NAV  T-beam  fuselage  fabrication  process:  (a)  25.4  mm  x25.4  mm 
xl  mm  finely  polished,  Cr/Au  electroded  PZT-5H  chip  (b)  machining  of  PZT  using 
high  precision  dicing  saw  (c)  patterning  of  bottom  electrode 


The  web  contracts  due  to  the  applied  field,  producing  a  differential  strain  across 
the  cross  section  and  constant  curvature  bending.  With  a  fixed-free  boundary 
condition,  a  cantilevered  T-beam  bends  out-of-plane  to  produce  relatively  large 
displacement. 

These  1  mm  thick,  20  mm  long  actuators  produce  200  /im  peak  to  peak  free  tip 
displacement  and  30  mN  blocking  force  at  1  V /  fim  DC.  T-Beams  are  fabricated 
using  a  K&S  980  Series  high  precision  dicing  saw  with  a  100  fim  wide  nickel  blade 
that  is  coated  in  3-6  fim  diamond  grit.  The  fuselage  is  diced  from  a  25.4  mm  x 
25.4  mm  x  1  mm  finely  polished,  Cr/Au  electroded  PZT-5H  chip  that  is  bonded 
to  a  2”  x  3”  glass  slide  using  1827  photoresist.  A  5  fim  layer  of  1827  spun  onto 
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Figure  2.2.2.  Schematic  diagram  of  the  PSNAV:(a)  Top  view,  (b)  Cross-Section  of  a 
T-Beam  Actuator 

the  PZT  chip  protects  the  T-Beams  from  particles  during  dicing.  The  fabrication 
process  steps  are  shown  in  section  2.2.  The  dicing  blade  spins  at  30  krpm  and 
machines  the  PZT  at  a  rate  of  0.3  mm/s.  Multiple  passes  realize  the  T  cross 
section,  and  create  a  20  mm  x  110  pm  air  gap  in  between  adjacent  T-beams. 
This  leaves  approximately  5  mm  of  PZT  material  at  the  base  to  ensure  that  the 
three  beams  remain  parallel  and  rigidly  connected.  The  center  beam  in  this  work 
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is  passive,  but  could  potentially  be  biased  and/or  sinusoidally  actuated  to  move 
the  thrust  vector  and  increase  flap  amplitude.  The  two  outer  T-Beams  create 
displacement  which  is  amplified  by  the  hinge  and  lever  system  producing  large 
angular  rotation.  The  T-beams  are  mounted  to  a  custom  PCB  board  using  silver 
epoxy,  and  wirebonded. 

The  schematic  model  of  the  NAV  system  is  shown  in  Figure  2.2.2.  Beams  1,  2 
and  3  are  the  piezoelectric  T-beam  actuators  and  beams  4  and  5  are  wing  supports. 
The  T-beams  are  fixed  at  one  end  and  connected  to  wing  supports  at  the  other 
end.  Beam  4  is  connected  to  beams  2  and  3  with  hinges  at  b-b’  and  d-d’  and 
beam  5  is  hinged  to  T-beams  2  and  3  at  d-d”  and  f-f”.  Cantilevered  T-beam 
models  developed  in  a  previous  paper  [60]  are  used  to  predict  the  displacement, 
u  and  the  blocking  force,  F.  The  parameters  that  define  the  cross-section  are  the 
total  height,  h,  total  width,  s,  web  thickness,  b,  and  flange  thickness,  t  as  shown 
in  Figure  2.2.2(b). 

T-beam  displacement  increases  with  the  square  of  the  length,  Lact  [60].  The 
current  design  has  constraints  on  the  length,  Lact  =  2  cm  and  total  height,  h  = 
1  mm  because  it  is  fabricated  from  a  25.4  mm  x  25.4  mm  x  1  mm  PZT  chip.  We 
constrain  s  =  1  mm  and  2  mm  for  the  outer  and  center  T-Beams,  respectively. 
The  center  T-Beam  is  twice  as  wide  because  it  is  loaded  by  both  outer  T-Beams 
and  requires  double  the  stiffness.  Web  thickness,  b,  and  flange  thickness,  t,  are 
chosen  to  optimize  the  T-beam  mechanical  energy  EMEch  =  uF.  Figure  2.2.3 
shows  the  maximum  mechanical  energy  output  for  a  given  web  thickness  and  its 
corresponding  optimal  flange  thickness  for  the  outer  T-beams.  The  optimal  flange 
and  web  thickness  are  331  pm  and  245  pm,  respectively.  Similarly,  the  optimized 
central  T-beam  has  t  =  331  fim  and,  b  =  490  pm. 
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Web  thickness,  b  (jum) 


Figure  2.2.3.  T-beam  maximum  mechanical  energy,  Emech  versus  web  thickness, 
6,  for  optimal  flange  thickness,  t,  (blue-solid)  and  T-beam  optimal  flange  thickness,  t, 
versus  web  thickness,  b,  (green-dashed)  length,  Lact  =  2  cm,  total  height,  h  =  1  mm, 
and  total  width,  s  =  1  mm 

2.2.2  Amplification  Mechanism 

Using  piezoelectric  actuators,  the  small  linear  displacement  must  be  amplified 
and  converted  into  large  angular  displacment  sufficient  for  flapping.  The  clap¬ 
ping  mechanism  outlined  in  this  work  amplifies  the  T-beam  displacement  using 
miniature  revolute  joints  and  hinges  spaced  a  certain  distance  apart,  creating  a 
lever  mechanism.  The  revolute  joints  are  fabricated  by  inserting  3-5  mm  long 
162.5  /mi  tungsten  pins  into  2  mm  long  165  /im  inside  diameter  medical  syringes 
housed  in  2  mm  silicon  extensions,  which  ensure  proper  alignment  and  bonding 
and  attach  to  the  free  end  of  the  T-beam  actuators,  shown  in  Figure  2.2.4. 


Optimal  flange  thickness,  t  (jum) 
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Syringe  Needle 


Tungsten  Rod  Silver  Epoxy 


Figure  2.2.4.  Drawing  of  revolute  joint  made  from  diced  1  mm  thick  Silicon,  165 
firn  inside  diameter  medical  grade  syringe  needle,  162.5  f.im  diameter  tungsten  rod,  and 
silver  epoxy  beads. 

2.2.3  Wing  Design 

The  wing  is  fabricated  by  first  stretching  12.5  /irn  thick  Mylar  sheet  onto  a  frame, 
then  precisely  placing  a  70  /im  thick  aluminum  foil  wingframe  onto  prescibed 
location  using  a  linear  stage.  The  wingframe  is  shaped  similar  to  a  fruitfly  aerofoil. 
The  wing  rod,  200  /im  x  600  [irn  x  6  cm  stainless  steel,  is  placed  in  a  similar  fashion 
and  connects  the  two  aerofoils  forming  the  leading  edge.  A  150  /im  gap  between 
the  leading  edge  and  the  wing  frame  forms  the  flexure  hinge  which  allows  passive 
wing  rotation.  Figure  2.2.5  shows  a  picture  of  a  fabricated  wing  set,  and  zoomed 
in  picture  of  the  flexure  hinge.  A  15°  bend  is  introduced  in  the  wing  rod  to  give 
clearance  for  the  second  wing  set.  The  wingsets  and  T-beams  are  set  on  separate 
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Figure  2.2.5.  100  / im  thick  aluminum  wing  frames  mounted  on  12.5  ^ m  thick  mylar 

wing  membrane  bonded  to  6  cm.  long  stainless  steel  wing  supports 
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Figure  2.2.6.  Schematic  of  fabricated  T-beam  fuselage,  revolute  joints,  and  wing  rod 
comprising  the  wing  stroke  amplification  mechanism 


micrometer  stages  and  the  leading  edge  of  the  wing  rod  is  precisely  connected  to 
the  tungsten  pins  under  a  microscope,  and  bonded  using  epoxy.  An  end  view 
of  the  fabricated  T-beams,  amplification  mechanism,  and  wing  rods  is  shown  in 
Figure  2.2.6. 
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2.3  Modeling 
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(c)  Front  View 


(d)  Right  Side  View 


Figure  2.3.1.  Schematic  of  the  clapping  winged  PSNAV  system  model  of  one  wing 
set. 

Figure  2.3.1  shows  a  schematic  drawing  of  the  clapping  wing  mechanism.  The 
goal  of  the  structural  dynamic  modeling  is  to  understand  the  resonance  behavior 
of  the  mechanism  to  aid  in  future  designs.  Thus  to  simplify  the  order  of  the 
model,  a  single  wingset  is  modeled,  and  the  interactions  with  the  other  wingset  is 
neglected.  The  mechanism  is  clamped  and  oriented  vertically  as  in  Figure  2.3.1(c). 
The  wing  rod  is  assumed  to  be  rigid  and  connects  the  leading  edge  of  each  aerofoil 
to  the  central  T-beam,  also  assumed  to  be  rigid,  at  the  hinge  at  the  point  O.  The 
flapping  stroke  is  constrained  to  the  xy  plane  and  approximated  for  small  angles 
by  sin  (0)  «  0  =  yx/La  where  yx  is  the  T-beam  displacement.  The  T-beam 
actuator  and  pin  assembly  is  modeled  as  a  force,  FT,  and  linear  spring,  kT,  hinged 
to  the  wing  rod  at  point  T.  The  aerofoils  A  and  B  attach  to  the  leading  edge  by 
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flexure  hinges  represented  as  hinged  joints  with  added  torsional  springs,  k^A  and 
kfe,  respectively.  The  leading  edges  of  aerofoils  A  and  B  undergo  displacement 
by  Rcg A and  respectively.  The  center  of  gravity  and  inertia  of  the  wings 

are  calculated  using  SolidWorks.  The  mass  of  the  T-beam  actuator  is  negligible 
compared  to  the  inertia  of  the  wings. 

The  sources  of  losses  in  the  design  are  in  the  revolute  joints  and  in  the  aero¬ 
dynamic  drag  experienced  by  the  aerofoils.  The  lossed  in  the  hinges  are  modeled 
as  a  single  damper,  applied  at  the  hinge  at  O.  To  properly  model  the  unsteady 
aerodynamics  and  the  fluid  structure  interation  is  beyond  the  scope  of  this  work. 
Alternatively,  a  linearized  aerodynamic  loading  effect  is  simulated  upon  aerofoils 
A  and  B  by  introducing  viscous  dampers,  c^A  and  c^B,  applied  at  the  hinges  A 
and  B ,  respectively. 

r  i T 

Choose  as  generalized  coordinates  q  —  ipA  <t>  i^B  >  where  ip  a  and  ips  are 
the  wing  rotation  angles  of  the  left  and  right  wing,  respectively,  and  <p  is  the  flap¬ 
ping  stroke  angle  as  shown  in  Figure  2.3.1.  Using  the  small  angle  approximation, 
the  kinetic  and  potential  energy  expressions  can  be  written  as 

T  =  t4>  +  =  iU/jq 

i=A,B 

v  =  u  +  E  v*  =  lqTKq  (2.i) 

i=A,B 

where 


J 


K 


Jtp,A  J(f>,ij),A  0 

J <f>  J(j>,iij,B 

0  J(f>,xl),B  Jip,B 

kipA  0  0 

0^0  , 

0  0  k^B 


(2.2) 


(2.3) 
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the  terms  J^,  for  i  e  { A ,  B }  are  calculated  from  Solidworks,  and 


kef,  —  krL2g , 

k^i  =  kfiexi  +  mF^gRcg^yf  or  i  e  { A ,  B} 


(2.4) 

(2.5) 


are  calculated  using  parameters  from  Table  2.3.1.  The  Solidworks  model  is  cali¬ 
brated  using  the  measured  mass  of  the  wing  set.  The  equations  of  motion  can  be 
written  using  the  Euler- Lagrange  formulation,  where  the  Lagrangian 


L  —  T  —  V  —  iqTJq  -  tqTKq 


and  Lagrange’s  Equations 


d  ( dL 
dt  \<9q 


where  the  generalized  forces  are 


dL 

~  -7T-  =  Jq  +  Kq  =  Q, 

oq 


Q  =  -Cq 


0 

FTLg 

0 


with  viscous  damping 


C  = 


0  C0  0 
0  0  ^ 


(2.6) 


(2.7) 


(2.8) 


and  FrLfJ  is  the  torque  generated  from  T-bearn  tip  force.  Substitution  of  the 
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generalized  forces,  intoEquation  2.7  produces  the  PSNAV  system  dynamics 


d 

dt 


q 

q 


o  i 
J  lK  j  xc 


q 

q 


J1 


o 

Q  AF 


(2.9) 


Table  2.3.1.  PSNAV  Parameters  and  Values 


Parameter 

Description 

Value 

R 

span  of  one  wing,  root  to  tip 

37.5 

mm 

c 

chord  of  wing 

17 

mm 

Lg 

distance  from  fulcrum  to  T-beam 

1.68 

mm 

L  flex 

flexure  length 

1.5 

pm 

Lact 

T-beam  length 

20.2 

mm 

Lpin 

Tungsten  pin  length 

1.63 

mm 

t mylar 

thickness  of  Mylar 

12.5 

pm 

Emylar 

Young’s  modulus  of  mylar 

4 

GPa 

EpZTbH 

Young’s  modulus  of  PZT5H 

63 

GPa 

It 

T-beam  area  moment  of  inertia 

0.0371 

mm 4 

I,,  _  GEpztshIt 

T-beam  stiffness 

817 

N 

1  E^ct(2Lact-\-Epin) 

m 

pAl 

density  of  Aluminum 

2700 

kg_ 
m 3 

Pss 

density  of  stainless  steel 

8000 

kg_ 
m 3 

Rcg,x 

aerofoil  CG  from  root,  spanwise 

17.5 

mm 

Rcg,y 

aerofoil  CG  from  root,  chordwise 

8.5 

mm 

mp 

wing  frame  mass 

116 

mg 

mR 

wing  rod  mass 

23 

mg 

T  ^  flex^w 

flex  ~  12 

flexure  area  moment  of  inertia 

3.3e-6 

mm 4 

P  Emyidj-I  fiex 

Eflex  r  ,7 

J  ^  flex 

flexure  stiffness 

0.088 

N  ■  mm 

Figure  2.3.2  shows  a  cross  section  of  the  wing  chord  illustrating  the  aerodynamic 
forces  and  torques  acting  upon  the  wing  and  the  definition  of  angle  of  attack, 
a  =  arctan  (ry/tv),  which  is  crucial  for  lift  production.  First,  the  translational 
velocity  of  the  wing  is  transformed  into  the  wing  reference  frame  by  the  rotation 
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Z 


Figure  2.3.2.  Aerodynamic  forces  on  the  wings(adapted  from  Andersen  et.  al.  [61]). 


matrix 


R(,0) 


—  sin  (-0)  cos  (ijj) 

—  cos  (j/j)  —  sin  (r/>) 


Following  the  approach  in  [61],  the  viscous  torque 


r 


V 


1 

64 


RpairC(n/2)c4\'tp\'ljj 


(2.10) 


and  forces 


P'x'  Pair^Vyf  Fxt 


(2.11) 
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Fy'  =  Pair ?VX>  -  Fy,  (2.12) 

act  on  one  aerofoil  in  the  wing  reference  frame  where  the  circulation  around  the 
aerofoil  is 

T  =  —-CTrRc\v\sin(2a)  +  - CrRc 2^, 
and  the  viscous  drag  force  is 

F1"  =  ^PairRc(CD( 0)cos2  (a)  +  Cd  ( 7J)  sin2  (a))|v|v'. 

The  drag  coefficient  Cp  is  evaluated  at  a  =  0  and  a  =  7t/2  using  data  from  [61]. 
The  torque  and  force,  transformed  back  into  the  inertial  reference  frame,  are 

Rc 

T  =  —Pair T|v|  -  t",  (2.13) 

F  =  R~1('0)F/,  (2.14) 

respectively.  The  inertial  forces  and  torques  in  Equation  2.13)  act  at  the  center  of 
pressure,  a  quarter  chord  from  the  leading  edge  of  the  aerofoil.  The  virtual  work 
associated  with  the  aerodynamic  loads  and  the  T-bearn  force,  FT,  is 

SW  =  FTSyT  +  ^2  {rMi  +  FijX8rijX  +  Fi)y5ri)y}  (2.15) 

i=A,B 

=  (FtLg+  ^2  Fi)XRcP}X)8(j)  +  ^2  {  (t*  -  Fiy|)  (2.16) 

i=A,B  i=A,B 

=  Q<t,S(j)+  ^2  Qpi^i  (2.17) 

i=A,B 

Substitution  of  the  generalized  aerodynamic  forces,  Q af  =  [QpAiQ<i»Qpb\T 
into  Eq.  ??  produces  the  PSNAV  system  dynamics 
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2.4  Experimental  Testing  and  Results 
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Figure  2.4.1.  Photograph  of  experimental  test  setup  for  (a)  T-beam  actuators  (b) 
nano  air  vehicle 

A  photograph  of  the  experimental  setup  for  actuator  measurement  is  shown  in 
Figure  2.4.1(a).  An  input  signal  is  generated  in  Lab  VIEW  and  amplified  using  a 
Trek  Model  609E-6  high  voltage  amplifier  which  can  provide  lOkV  max  voltage 
up  to  6kHz.  A  Polytec  Laser  Vibrometer®  with  a  10X  objective  focuses  a  20  /im 
laser  onto  the  gold  web  electrode  of  a  T-Bearn  and  returns  displacement  signals  to 
Lab  VIEW.  A  The  clapping  wing  mechanism  is  mounted  vertically,  and  Newport 
XYZ  linear  stages  allow  precision  positioning.  The  wing  trajectories  of  the  PSNAV 
are  measured  using  strobed  photography.  The  PSNAV  is  fixed  to  an  ABS  plastic 
lift  apparatus  that  uses  a  razor  blade  fulcrum  that  preloads  and  transmits  dynamic 
loads  to  an  Aurora  Scientific,  Inc.  402A  force  transducer  with  10  /iN  resolution 
and  2kHz  bandwidth.  At  least  50  cycles  are  measured  to  obtain  steady-state  lift. 

The  T-Beam  actuators  are  tested  off-resonance,  quasi-statically  at  1  Hz. Figure 


28 


Figure  2.4.2.  Quasi-static  free  tip  displacement  of  a  T-beam  actuator 

2.4.2  shows  the  free  tip  displacement  and  blocking  force  response  of  the  left  and 
right  T-Beam  under  varying  peak  applied  fields.  The  actuators  are  designed  to 
provide  large  displacement,  so  the  hysteresis  associated  with  relatively  high  strain 
PZT-5H  is  a  minor  drawback.  During  blocking  force  measurements  the  compliance 
of  the  PCB  board  was  apparent.  The  measured  blocking  force  is  significantly  less 
than  the  true  blocking  force  of  the  T-beam  actuator  due  to  bending  of  the  PCB 
board  seen  in  figure  2.4.3.  Ideally,  the  displacement  at  the  base  of  the  T-beam 
would  be  zero  during  blocked  force  measurement. 
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Figure  2.4.3.  Measured  force  of  T-beam  actuator.  Actual  blocking  force  is  larger  due 
to  compliance  in  mounting. 

A  field  of  0.7  V / [im  amplitude  with  0.1  V / jim  DC  bias  is  applied  to  the  fabri¬ 
cated  prototype  PSNAV  from  DC  to  30  Hz.  Two  resonance  modes  were  observed 
near  9.5  Hz  and  25.5  Hz.  Figure  2.4.4  shows  a  photograph  of  the  PSNAV  flapping 
stroke  motion  in  the  open  and  closed  position,  showing  ~11  —  14°  peak  to  peak 
per  aerofoil  at  9.5  Hz  that  produces  near  clapping  (Figure  2.4.4(a))  and  a  gap  of 
around  30°  (Figure  2.4.4(b))  at  different  times  during  the  flapping  cycle.  Figure 
2.4.5  shows  images  at  five  stages  of  the  flapping  cycle  at  9.5  Hz.  The  left  set  of 
wings  passively  rotate  from  -33°  to  21°  during  the  flapping  cycle.  The  right  set 
of  wings,  however,  show  substantially  less  rotation  and  translation.  This  is  due 
to  the  slight  differences  in  resonance  frequencies  between  the  two  sets  of  wings. 
Figure  2.4  show  the  flapping  and  rotation  amplitudes  versus  frequency.  The  right 
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Figure  2.4.4.  LionFly  in  the  closed(a)  and  open(b)  position  showing  30  deg.  of  stroke 
from  DC  to  9.5Hz  using  0.7  V/gm  amplitude  and  0.1  V/gm  DC  bias 


Figure  2.4.5.  Strobed  images  of  one  period  of  wing  rotation  actuated  at  9Hz  using 
0.7  V/gm  amplitude  and  0.1  V/gm  DC  bias 

pair  of  wings  are  slightly  off-resonance  so  the  amplitudes  are  smaller.  The  reso¬ 
nance  mode  at  25.5  Hz  produces  ~8°  peak  to  peak  flapping  angle  and  ~24°  rotation 
angle.  The  response  of  4>a,b  is  nearly  90°  out  of  phase  with  6  at  9.5  Hz  and  180°  at 
25.5  Hz.  Figure  2.4.7  shows  the  experimentally  measured  thrust  versus  frequency 
at  0.1  V/um  bias  and  0.7  V/um  amplitude.  The  peaks  at  9.5  Hz  and  25.5  Hz  are 
clearly  seen. 
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Figure  2.4.6.  Stroboscopically  photographed,  image  processed  experimental  data  and 
simulated  wing  A  rotation  (top),  flapping  </>  (middle),  and  wing  B  rotation,  i^b  (bot¬ 
tom)  versus  frequency  at  applied  field  of  0.7  V I [im  and  0.1  V / fim  DC  bias:  experimental 
data  for  wing  set  1  (blue  circles),  simulated  data  for  wing  set  1  (blue  solid),  experimental 
data  for  wing  set  2  (green  squares),  and  simulated  data  for  wing  set  2  (green  dashed). 
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Figure  2.4.7.  Average  lift  force  of  NAV  versus  frequency  under  applied  field  of  0.7 
V / firri  and  0.1  V / /im  DC  bias.  Resonance  frequencies  of  the  two  modes  are  highlighted 
(red  dashed). 


Chapter  3 

Monolithic  SUEX  Flapping  Wing 
Mechanisms  for  Pico  Air  Vehicle 
Applications 

3.1  Introduction 

This  chapter  introduces  a  simple  flapping  wing  mechanism  fabricated  monolithi- 
cally  from  SUEX  dry  film,  an  epoxy  based  negative  photoresist  similar  to  SU-8.  We 
present  the  fabrication  process  and  characterization  of  SUEX,  static  and  dynamic 
modeling  of  the  PAV,  and  experimental  results  from  multiple  prototype  devices 
that  demonstrate  applicability  of  the  developed  fabrication  process  for  flapping 
wing  air  vehicle  applications. 
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Figure  3.1.1.  Conceptual  drawing  of  the  Penn  State  LionFly:  (a)  Isometric  view, 
(b)  zoomed  inset  of  compliant  mechanism  and  piezoelectric  bimorph  actuator,  and  (c) 
operation  of  compliant  flapping  wing  mechanism  with  actuator  displacement  input,  x , 
and  flapping  output,  q b. 


3.2  Material  Choice  and  Fabrication 

3.2.1  Material  Choice  for  Compliant  Mechanisms 

Compliant  mechanisms  typically  consist  of  a  system  of  rigid  links  connected  by 
flexible  members  designed  to  bend  within  a  certain  range  of  motion.  These  mech¬ 
anisms  can  be  fabricated  from  multiple  layers  bonded  together  or  monolithically 
from  a  single  layer  with  an  etched  flexible  region.  The  use  of  multiple  layers 
can  increase  fabrication  complexity  and  introduce  problems  such  as  delamination. 
Monolithic  fabrication  is  simpler,  but  it  is  important  to  select  a  material  that 
can  be  fabricated  to  the  desired  geometry  with  suitable  material  properties.  Fig¬ 
ure  3.2.1  shows  a  typical  flexure  hinge  under  normal  operation,  where  one  rigid 
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Figure  3.2.1.  Schematic  of  a  flexure  hinge:  a)  dimensions,  b)  typical  operation 

link  is  fixed,  and  the  other  is  free  and  has  a  force  perpendicular  to  the  tip  of  the 
flexure  region.  The  bending  stiffness  of  the  rigid  portion  is  much  larger  than  the 
bending  stiffness  of  the  flexure  region  ( EIl  »  Elf).  This  is  achieved  by  reducing 
the  thickness  of  the  flexure  region,  in  which  a  reduction  by  a  factor  of  5  creates  a 
bending  stiffness  ratio  of  125  for  a  monolithic  flexure  hinge  because  the  moment  of 
inertia  for  a  rectangular  cross  section,  /  =  wt3  / 12.  The  maximum  bending  angle 


36 


for  a  short  rotational  flexure  can  be  determined  by  [19] 

=  2-^rL  PU 

YU/ 

where  ay  is  the  yield  strength  and  Ef  is  the  Young’s  modulus  of  the  material.  The 
factor  cry/Ef  determines  the  maximum  strain  before  yielding,  ey,  of  a  material,  in¬ 
dependent  of  geometry.  Table  3.2.1  shows  many  candidate  flexure  hinge  materials 
and  their  properties  ranked  by  yield  strain.  Ranking  high  are  polymer  materials 
such  as  polyimide,  and  low  are  metals  such  as  steel  and  polysilicon.  Single  crystal 
silicon  is  an  excellent  material  for  compliant  mechanisms  in  theory.  However  in 
practice,  sensitivity  to  crystallographic  defects  and  surface  imperfections  make  sin¬ 
gle  crystal  silicon  susceptable  to  crack  propagation  [62]  and  incapable  of  handling 
large  loads  [17]. 

SU-8,  a  negative  tone,  epoxy  based,  near  UV,  thick  photoresist,  has  tradition¬ 
ally  been  used  for  high  aspect  ratio  applications  such  as  electroplating  masks  [63] , 
deep  reactive  ion  etching  masks  [64],  and  microfluidics  [65].  Its  use  as  a  structural 
material  has  been  shown  in  applications  such  as  compliant  microgrippers  [59]  and 
recently  in  NAV  applications  [21,  10].  SU-8  is  typically  spin  coated  onto  a  Si  wafer, 
soft-baked  (65-95  °C),  exposed  to  365  nm  UV  light,  post-exposure  baked  (PEB) 
(65-95  °C),  developed  in  PGMEA,  and  optionally  hard  baked  (HB)  (150-200  °C). 
Multiple  layers  must  be  spun  after  the  initial  layer  soft  bake  step  to  increase  over¬ 
all  thickness.  SUEX,  by  DJ  DevCorp,  Inc.  has  nearly  identical  characteristics  to 
SU-8  [66],  but  is  available  as  a  dry  him  in  thicknesses  up  to  1  mm.  Post  soft  bake, 
the  processing  is  nearly  identical  to  SU-8,  but  SUEX  can  be  processed  as  a  sub¬ 
strate,  with  the  exception  of  50  /im  PET  protective  layers  that  are  easily  removed. 
Low  density  and  fabrication  advantages  make  SUEX  an  excellent  material  for  PAV 
applications,  despite  not  ranking  first  in  yield  strain  in  Table  3.2.1. 
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Table  3.2.1.  Table  of  candidate  materials  for  flexure  Hinges 


Material 

ey  (%) 

E  ( GPa ) 

<7y  (MPa) 

p  ( kg/m 3) 

Source 

PI-2525 

5.35 

2.4 

128.5 

1420 

[67] 

Si“ 

3.68 

190.0 

7000.0 

2300 

[62] 

Si3N4“ 

3.64 

385.0 

14000.0 

3100 

[62] 

SUEX 

3.07 

2.8 

86.0 

1190 

[68] 

PET 

2.87 

2.8 

79.3 

1356 

[69] 

PP 

2.43 

1.4 

34.0 

855 

[70] 

Parylene  C 

1.99 

2.8 

55.0 

1289 

[71] 

Nylon-66 

1.96 

2.8 

55.0 

1150 

[70] 

Ti-136 

1.03 

114.0 

1170.0 

4820 

[70] 

Steel  4140c 

0.79 

207.0 

1641.0 

7750 

[70] 

A1  70756 

0.70 

71.7 

503.0 

2795 

[70] 

poly-Si 

0.55 

169.0 

930.0 

2329 

[70] 

“single  crystal,  6heat  treated,  CQ&T400 


3.2.2  Monolithic  Multilevel  Fabrication  of  SUEX 

Similar  to  other  negative  photoresists  [72,  73],  SUEX  has  unique  optical  properties 
where  absorption  is  much  higher  at  shorter  wavelengths,  allowing  only  partial 
penetration  of  thickness.  At  higher  wavelengths  (above  365  nm)  the  material  is 
transparent,  allowing  ultra  thick  and  high  aspect  ratio  features.  Figure  3.2.2  shows 
the  optical  transmission  spectrum  of  unexposed  SUEX  film,  taken  using  a  Perkin 
Elmer  Lambda  950  UV/Vis  Spectrometer.  The  film  is  nearly  transparent  above 
365  nm,  but  blocks  wavelengths  below  350  nm,  enabling  definition  of  thin  flexible, 
full  thickness,  and  through  etched  regions.  We  exploit  this  phenomenon  to  define 
flexure  hinges  using  UV  energy  at  310  nm  for  thin  flexible  regions,  and  365  nm 
for  thick  rigid  regions.  Figure  3.2.3  shows  the  fabrication  process:  (a)  250  prri 
thick  SUEX  comes  sandwiched  between  two  50  micron  PET  layers  for  support.  In 
(b)  the  top  layer  is  removed,  and  the  first  mask  defines  the  thin  flexure  and  wing 
membrane  region  with  the  thickness  controlled  by  the  exposure  dose  at  310  nm 
wavelength.  Figure  3.2.3  (c)  shows  the  thick  or  rigid  link  regions  of  the  LionFly 
defined  by  exposure  at  365  nm  using  the  second  mask.  The  post  exposure  bake 
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(PEB)  (d),  performed  at  65°C  for  30  minutes  on  a  hotplate,  is  crucial  for  cross 
linking  of  SUEX.  The  wafer  is  removed  and  allowed  to  cool  for  an  additional  half 
hour.  The  bottom  PET  layer  is  removed,  and  the  SUEX  is  developed  in  a  mesh 
cage  in  PGMEA  (e)  for  one  hour  at  room  temperature  with  light  agitation.  The 
devices  are  removed  from  the  mesh,  blow  dried  with  jV2,  and  put  in  a  vacuum  oven 
at  50  °C  for  a  half  hour  to  dry.  The  devices  at  this  point  are  ready  for  mounting 
with  total  clean  room  processing  time  approximately  three  hours. 

Multiple  tests  at  varying  doses  were  conducted  to  characterize  the  effect  of 
310  nrri  exposure  dose  on  flexure  thickness.  The  UV  exposure  wavelength  can  be 
controlled  by  installing  different  filter  sets  in  mask  aligners  with  broadband  light 
intensity.  In  this  work,  two  mask  aligners  with  different  UV  intensity  at  310,  365 
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a)  250  pm  SUEX,  50  pm  PET  d)  Post  Exposure  Bake 


b)  Mask  1:310  nm  Exposure  e)  PGMEA  Development 


Through  Etch  Flexure  Rigid 

c)  Mask  2:  365  nm  Exposure 

PET  Substrate 
Unexposed  SUEX 
Exposed  SUEX 
Crosslinked  SUEX 


Figure  3.2.3.  Fabrication  of  SUEX  structure:  (a)  initial  SUEX  substrate  between  two 
50  fim  PET  layers,  (b)  310  nm  exposure  dose  of  500  mJ/cm 2  for  25  firri  flexure  layer, 
(c)  365  nm  exposure  dose  of  1350  mJ/cm2  for  250  pm  rigid  layer,  (d)  post  exposure 
bake  (65°for  30  min)  to  activate  cross  linking  in  SUEX,  (e)  development  in  PGMEA  for 
1  hr  with  light  agitation  releases  structures 


and  405  nm  wavelength  were  used.  The  Karl  Suss  MJB-3  contact  aligner  was  used 
for  flexure  definition  and  has  average  light  intensity  of  2.9  mW/cm2  at  310  nm, 
near  zero  at  365  nm,  and  1.15  mW/cm 2  at  405  nm.  The  EVG  620  mask  aligner 
was  used  for  rigid  link  definition  and  has  average  light  intensity  of  4.5  mW/cm 2  at 
310  nm,  13.96  mW/cm2  at  365  nm,  and  26.14  mW/cm 2  at  405  nm.  Step  height 
measurement  was  conducted  using  a  Zygo®  NewView  100C  optical  prohlometer, 
and  total  thickness  measurements  using  an  optical  microscope.  Figure  3.2.4  shows 
the  flexure  thickness  begins  to  saturate  at  approximately  70±5  nm  at  doses  beyond 
1  J /cm2 ,  and  is  nearly  linear  between  0.25-0.75  J /cm2.  Figure  3.2.5  shows  a  cross 
section  of  a  fabricated  flexure  hinge  from  SUEX.  It  should  be  noted  that  in  this 
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Figure  3.2.4.  Flexure  thickness  versus  310  nm  exposure  dose 

work  only  two  thickness  levels  are  chosen,  however  this  process  can  be  extended  to 
have  multiple  levels  with  additional  masks  defining  each  level  with  increasing  310 
nm  exposure  dose. 

3.2.3  Actuator  Selection 

In  this  work,  we  use  high  performance  PZT-5H  bimorphs  from  Piezo  Systems, 
Inc.®,  comprised  of  a  115  gm  thick  composite  shim  material  sandwiched  between 
two  127  fim  thick  PZT-5H  piezoelectric  layers.  The  material  is  diced  into  smaller 
2  mm  x  10  mm  bimorphs  using  a  K&S  Model  980  dicing  saw.  Typical  actuators 
weigh  55  mg,  produce  ±  100 gm  tip  displacement  and  45  mN  blocking  force  at 
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Figure  3.2.5.  Cross  section  of  a  fabricated  flexure  hinge 


±150  V,  and  have  a  resonance  frequency  of  >  1  kHz. 

3.3  Compliant  Mechanism  Design  and  Modeling 

3.3.1  Mechanism  Design 

Due  to  the  flexure  length  limitations  set  by  Equation  3.1,  it  is  difficult  to  achieve 
large  amplification  with  out-of-plane  displacement.  In  this  work,  an  offset  slider- 
rocker  compliant  mechanism  is  designed  using  a  bimorph  mounted  such  that  the 
input  displacement  direction  is  in-plane.  This  mechanism  has  high  amplification 
and  is  shifted  from  the  horizontal  plane  using  a  small  vertical  offset.  To  realize 
this  design,  flexure  hinges  must  be  designed  to  undergo  large  deformation  yet  resist 
buckling  when  loaded.  Using  the  standard  Euler  buckling  relationship  for  a  flexure 
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Figure  3.3.1.  Schematic  of  pseudo  rigid  body  model  of  slider  crank  mechanism 


that  is  fixed  at  one  end  and  free  on  the  other,  the  critical  axial  load  is 


Fcrit 


i f-Ejwjt) 

48  L) 


(3.2) 


In  this  work,  the  flexure  lengths  are  designed  to  bend  at  a  maximum  of  60°  and 
the  flexure  widths  are  designed  to  prevent  buckling. 


3.3.2  Modeling 

Figure  3.3.1  shows  a  schematic  of  the  psuedo  rigid  body  approximation  of  the 
compliant  slider-rocker  mechanism.  This  approximation  has  been  developed  for 
large  range  of  motion,  and  allows  the  flexure  hinges  to  be  approximated  as  ideal 
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pins  with  torsional  stiffnesses  of  [70] 


k 


eq 


lKe 


ESwft f 

12  Lf  ’ 


(3.3) 


where  7  is  the  characteristic  radius  factor  and  Kq  is  a  stiffness  coefficient  dependent 
on  the  range  of  motion.  This  equation  is  used  when  the  flexure  length  is  on  the 
same  order  of  magnitude  as  the  length  of  the  rigid  links.  It  is  also  assumed  that 
no  off  axis  translation  or  rotation  occurs  and  the  stiffness  is  linear  along  the  range 
of  motion  of  the  spring.  The  mechanism  consists  of  three  main  hinges,  O,  A,  and 
B ,  with  torsional  stiffness,  ko,  kA,  and  ks  respectively.  Linear  viscous  damping  is 
included  on  main  wing  hinge  O, 


Co  Coj  T  Co,air 


(3.4) 


where  Coj  is  a  structural  damping  term,  and  Co,air  is  linear  damping  term  to 
simulate  aerodynamic  drag,  with  Co,air  —  0  in  vacuum.  The  bimorph  is  modeled 
as  a  linear  voltage  controlled  force  acting  upon  the  slider  [74], 


Fact  —  GnrtV 


act  ' 


9EIeffY(2Y  +  l)d31V 
2  tm(tm  +  2tp)Lact(AlA3  +  3  Y2  +  3  Y  +  1) 


(3.5) 


with  linear  stiffness,  kact,  viscous  damping,  cact,  and  linear  displacement  a;. The 
actuator  stiffness  is 


kact 


3  EL 


eff 


T3 

^  act 


(3.6) 


with 


EIeff  = 


LCact  Ep 


r+3 


t  F  i3 

A  +  +  (t  +t?  +  —  — 
3  +  tp[Tm  +  p)  +  E„  6 


2 


(3.7) 
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Figure  3.3.2.  Link  diagram  of  the  slider  rocker  mechanism  (initial  (grey  dashed)  and 
final  (black  solid)  configurations) 

X  =  Em/Ep,  and  Y  =  tm/(2tp).  The  bimorph  damping  is  included  for  generality 
but  can  be  safely  neglected  relative  to  the  large  wing  damping.  The  bimorph  is 
attached  to  the  wing  by  the  coupler  link  with  stiffness  kB-  Figure  3.3.2  describes 
the  flapping  angle,  q 1,  defined  with  respect  to  the  negative  rc-axis,  the  hinge  B  angle, 
0B,  defined  with  respect  to  the  positive  x-axis,  the  relative  angle  between  the  links 
AO  and  AB,  6a,  and  the  actuator  input,  x.  The  psuedo  rigid  body  approximation 
gives  r2  =  Loa  +  ^Lo  ~  (1  —  7 )La,  where  Lo  and  La  are  the  lengths  of  the  flexure 
hinges  O  and  A  respectively.  The  loop  closure  equations  are  solved  to  yield 


r3((j),x) 

r2  cos (j)  —  ri(x) 
cos  6  b 

(3,8) 

9B{(t),x) 

ro  sin  (j)  —  r 4 

=  arctan - - — 

r2  cos  0  —  r  1  (x) 

(3.9) 

9a{4>,x) 

=  6b  -  0 

(3.10) 

where  ri(x)  =  r\0  +  x,  and  r±0  is  the  initial  length  of  r\(x). 
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Using  Figure  3.3.1,  the  potential  energy  of  the  mechanism  is 

U  =  ^0(0  -  4>ref)2  +\  X )  -  ^*,re/)2  (3.11) 

i=A,B 

^ (^3 (0;  *^)  ^3, re/)  ~h  ~^kact(x  Xre/)  , 

where  0re/,  9 a, ref,  9b, ref,  r3,ref ,  and  xre/  are  the  reference  positions  of  0,  0,4,  0b, 
r3,  and  x,  respectively.  Gravitational  potential  energy  is  neglegible  compared  to 
the  strain  potential  energy  stored  in  the  hinges  and  actuator. 

The  kinetic  energy  is 

T  =  ^ Jw4 J  +  ^rneqx2,  (3.12) 

where  Jw  is  the  wing  inertia  and  meq  is  the  equivalent  mass  of  the  actuator, 
0.23mact.  The  mass  and  inertia  of  the  coupler  link  is  neglibly  small.  The  wing 
inertia  is  approximated  from  a  computer  aided  design  (SolidWorks®)  model  of 
the  wing  with  the  density  adjusted  to  match  the  measured  mass  of  the  fabricated 
wing.  Assuming  the  wing  is  thin  flat  plate,  the  added  mass  can  be  found  by 


dm 


am 


7T 


jPo.i  rc(r)2dr, 


(3.13) 


where  pa%r  is  the  density  of  air,  and  c(r)  is  the  chord  length  at  position  r  along 
the  wing.  The  added  mass  is  calculated  by  numerically  integrating  this  quantity 
in  Matlab  using  a  discretized  image  of  the  wing.  Similarly,  the  added  inertia  is 
obtained  by  numerically  integration  of 

Jam.  =  /  r2 dm am  =  - pair  /  c(r)2r2dr.  (3.14) 

00  4  00 


The  added  mass  is  zero  in  vacuum. 
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The  virtual  work  associated  with  the  bimorph  actuator  and  damping  is 

5W  =  —c0<j>S(j)  +  (■ GactV  -  cactx)5x.  (3.15) 

Choosing  generalized  coordinates  q  =  [0,  x]T ,  the  equations  of  motion  can  be 
derived  using  the  Euler- Lagrange  formulation  as 


Jq  +  Cq  +  f(q)  =  Q, 


where 


J  = 

C  = 


Jw  0 
0  meq 

c0  0 
0  cac.t 


/(<?) 


9U_ 

d<p 


dx 


Q  = 


o 

GactV 


The  static  equilibrium  point,  qeq  depends  on  voltage  as  follows: 


fig)  -Q  =  o. 


(3.16) 


(3.17) 

(3.18) 

(3.19) 

(3.20) 


(3.21) 


The  nonlinear  stiffness  can  be  linearized 


dfi 

dfi 

d<t> 

dx 

d/2 

d/2 

d<f> 

dx 

<7=<Jeq 


kn  k\2 

k\2  k22 


(3.22) 


Taking  the  Laplace  transform  of  the  linearized  equations  of  motion  produces  the 
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Mechanism 


Flexure  Hinges 


PZT  Actuator 


r 

1  mm 


5  mm 


Figure  3.4.1.  Photograph  of  LionFly  prototype  07  and  zoomed  inset  of  flexure  hinges 
and  bimorph  connection  to  coupler  link. 


linearized  system  transfer  functions 

_ (  ki2Gact/ (^JwTTleq)') 

W)=—  DW 

X (s)  (5  (cq /  Jyj)s  “f"  k\\/ Jw)Gact/ ITleq) 

V(s)=  D(T) 


(3.23) 

(3.24) 


where 


D(s)  —  (s2  +  ( Co / Jw)s  +  kn/  Jw) 

(s  ~t~  (cacf/,m.eg)s  4“  k'22 / TTleq)  ky^/ wTTleqj . 


(3.25) 
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3.4  Experimental  Validation 

The  single  wing  LionFly  prototype  LF07,  shown  in  Figure  3.4.1,  is  used  for  model 
validation.  The  piezoelectric  bimorph  tip  is  bonded  to  the  coupler  link,  shown  in 
the  inset  of  Figure  3.4.1  using  precision  linear  stages.  The  actuator  is  clamped  at 
the  base  using  a  2  mm  by  2  mm  glass  square.  Lead  wires  are  soldered  directly 
onto  the  bimorph.  Gold  film  reflective  test  surfaces  are  attached  to  the  wing  and 
bimorph  for  laser  vibrometer  measurement. 

3.4.1  Experimental  Setup 

Figure  3.4.2  shows  the  experimental  setup  used  to  test  the  LionFly  mechanism. 
The  mechanism  is  mounted  on  a  carbon  fiber  holder  in  an  aluminum  vacuum  cham¬ 
ber  with  an  transparent  quartz  cover.  The  vacuum  chamber  is  fixed  to  Newport 
linear  stages  for  precise  alignment.  A  Polytec®  laser  vibrometer  measures  the 
displacement  at  the  test  points  on  the  wing  surface  and  the  bimorph  tip.  The 
vibrometer  uses  a  10  x  objective  lens  to  focus  the  laser  to  a  25  [mi  point  and  re¬ 
turns  displacement  and  velocity  data  to  LabView®  with  sub  micron  precision.  For 
small  signal  actuation,  the  LabView®  DAQ  system  can  provide  a  10V  amplitude 
signal,  but  for  high  voltage  experiments  (150  V  amplitude  with  150  V  DC  bias), 
a  Matsusada®  AP-1B3  amplifies  the  DAQ  output. 

3.4.2  Experimental  Results 

The  static  response  is  measured  using  an  applied  100  V  triangle  wave  over  240 
seconds.  The  theoretical  results  are  found  by  numerically  solving  the  nonlinear 
equilibrium  Equation  3.21  using  Newton’s  method.  The  theoretical  static  response 
in  Figure  3.4.3  (blue  solid)  matches  the  experimental  response  (blue  circles)  for 
small  voltages  but  is  lower  than  the  experimental  response  at  high  voltages  because 
PZT-5H  is  nonlinear,  showing  a  field  dependent  d^i  and  hysteresis  [75].  Introducing 
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Figure  3.4.2.  Photograph  of  experimental  testing  setup. 


a  quadratic,  field  dependent 

d31(V/t)  =  d31(V/t )  +  aiV/tf  (3.26) 

coefficient  produces  a  more  accurate  theoretical  (green  dashed)  curve.  Static  op¬ 
eration  shows  kinematic  amplification  of  0.244  deg./fim  and  low  voltage  gains  of 
7.3  deg./kV  and  31.24  /. im/kV  for  6  and  x  respectively. 

The  flap  and  actuator  displacement  frequency  responses  of  the  LionFly  mecha¬ 
nism  LF07  are  shown  in  Figures  3.4.4  and  3.4.5,  respectively.  Testing  results  in  air 
and  in  a  vacuum  chamber  pumped  down  to  76  Torr  show  velocity  data  measured 
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Voltage  ( V ) 


Figure  3.4.3.  Quasi-static  response  of  LF07:  experimental  data  (blue  circles)  ,  theoret¬ 
ical  curve  using  linear  d^i  (blue  solid),  theoretical  curve  using  quadratic  (green  dashed). 

from  multiple  points  on  the  wing.  The  parameters  used  in  model  validation  are 
shown  in  Table  3.4.1  for  the  bimorph  and  Table  3.4.2  for  the  LionFly  mechanism 
LF07.  Figure  3.4.4  shows  the  flapping  angle  calculated  from  the  outermost  test 
point  using  </>  =  arctan (zmeas/Rtp).  The  effect  of  aerodynamic  added  mass  and 
drag  are  seen  here  as  a  shift  in  resonance  frequency  and  a  reduction  in  amplitude, 
respectively.  The  theoretical  results  show  that  the  added  mass  calculation  is  ef¬ 
fective  in  capturing  this  effect.  The  linear  aerodynamic  damping  term  is  sufficient 
to  capture  the  reduction  in  amplitude  at  this  relatively  low  voltage  excitation.  No 
significant  wing  deformation  is  observed  using  the  data  from  the  numerous  test 
points  on  the  wing. 
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Figure  3.4.4.  Frequency  response  4>(s)/F(s)  of  LF07:  Experimental  data  in  air  (blue 
circles),  theoretical  response  in  air  (blue  solid),  experimental  data  in  vacuum  (green 
diamonds),  theoretical  response  in  vacuum  (green  dashed). 

The  frequency  response  of  the  LF07  bimorph  actuator  displacement  shows  good 
agreement  between  the  experimental  and  theoretical  results.  The  presence  of  the 
zero  in  Equation  3.24  is  clearly  seen  here,  at  53  Hz ,  and  shows  the  effect  of  the 
coupler  link  compliance. 

Several  LionFly  prototypes  were  tested  with  high  voltage  at  resonance,  and 
a  summary  of  their  performance  is  listed  in  Table  3.4.3.  The  fabricated  LionFly 
prototypes  show  a  consistent  resonance  frequency  of  49.13  ±  1.43  Hz  over  four 
devices.  Figure  3.4.6  shows  the  55°  peak  flapping  angle  of  the  LionFly  LF09.  Con¬ 
tinuous  operation  at  large  flapping  angles,  however,  led  to  failure  of  the  wing  hinge. 
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Figure  3.4.5.  Frequency  response  X(s)/V(s)  of  LF07:  Experimental  data  in  air  (blue 
circles),  theoretical  response  in  air  (blue  solid),  theoretical  response  in  vacuum  (green 
dashed). 


Longer  cycle  life  can  be  achieved  by  increasing  flexure  length  and/or  reducing  flex¬ 
ure  thickness  to  increase  the  maximum  bending  angle. 
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Table  3.4.1.  Bimorph  parameters  and  values  for  LF07 


Parameter 

Description 

Value 

dz  i 

piezoelectric  constant 

-274.0 

pm/V 

a 

quadratic  field  dependent  coefficient 

450.0 

( pm/mV )2 

EP 

Young’s  modulus  of  PZT-5H 

63.0 

GPa 

Em 

Young’s  modulus  of  shim 

2.0 

GPa 

Pp 

density  of  PZT-5H 

7500.0 

kg / m 3 

Pm 

density  of  shim 

1600.0 

kg /m3 

Lact 

free  length  of  actuator 

9.5 

mm 

tp 

thickness  of  PZT-5H  layer 

135.0 

pm 

t"m 

thickness  of  shim 

115.0 

pm 

W act 

width  of  actuator 

2.0 

mm 

TTlact 

mass  of  actuator 

54.0 

mg 

Figure  3.4.6.  Flapping  angle  of  LF09  with  150  V  amplitude,  150  V  DC  offset,  at  51 
Hz 
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Table  3.4.2.  LionFly  parameters  and  values  for  LF07 


Parameter 

Description 

Value 

Constants 

Ef 

Young’s  modulus  of  SUEX 

2.80 

GPa 

Pair 

density  of  air 

1.28 

kg/m 3 

PSUEX 

density  of  air 

1190 

kg/m 3 

7 

characteristic  radius  factor 

0.8517 

Ke 

stiffness  coefficient 

2.6762 

Mechanism  Parameters 

Lab 

length  of  coupler  link 

3.25 

mm 

Loa 

length  of  rocker  link 

520.00 

pm 

La 

length  of  flexure  hinge  A 

670.00 

pm 

Lb 

length  of  flexure  hinge  B 

275.00 

pm 

Lo 

length  of  flexure  hinge  0 

670.00 

pm 

tA 

thickness  of  flexure  hinge  A 

32.43 

pm 

tB 

thickness  of  flexure  hinge  B 

34.41 

pm 

to 

thickness  of  flexure  hinge  0 

33.02 

pm 

WA 

width  of  flexure  hinge  A 

2.30 

mm 

WB 

width  of  flexure  hinge  B 

2.30 

mm 

Wo 

width  of  flexure  hinge  0 

5.60 

mm 

( firef 

reference  position  of  (j) 

0.00 

rad. 

d  A, ref 

reference  position  of  9a 

0.00 

rad. 

@B,ref 

reference  position  of  9b 

7 r 

rad. 

Xref 

reference  position  of  x 

98 

pm 

rlo 

initial  value  of  r\ 

2.87 

mm 

r2 

rocker  length 

0.99 

mm 

Dj,  ref 

reference  value  of  r:i 

3.86 

mm 

U 

vertical  offset  of  slider 

0.10 

mm 

kL 

coupler  link  stiffness 

17.50 

N/mm 

Wing  Parameters 

mw 

mass  of  wing 

18.70 

mg 

Rw 

distance  from  hinge  0  to  wing  tip 

23.00 

mm 

RtP 

distance  from  hinge  0  to  test  point 

16.25 

mm 

J w 

wing  inertia,  in  vacuum 

1.80 

g  ■  mm2 

J am 

inertia  due  to  added  mass 

0.12 

g  ■  mm 2 

°OJ 

structural  damping  of  0 

8.80 

mN  ■  pm  ■  s 

C(D, air 

aerodynamic  damping  of  </> 

7.00 

mN  ■  pm  ■  s 

55 


Table  3.4.3.  Summary  of  LionFly  prototype  performance  at  150  V .  150  Vdc 


Device 

0(0)  (deg.) 

fr  (Hz) 

0Mai  (deg.) 

LF07 

0.16b 

50. 0b 

42c 

LF08 

/ 

47.5“ 

/ 

LF09 

12.1 

51.0 

55.0 

LF12 

10.5 

48.0 

30.0 

“100  Vpp,  device  failure  before  0  response  measured 
b20  Vpp ,  cIn  vacuum,  250  Vpp ,  125  V^c,  44  Hz 


Chapter  4 

Wing  Rotation  in  SUEX  Flapping 
Wing  Mechanisms 

4.1  Introduction 

In  this  chapter,  a  flexure  hinge  along  the  span  of  the  LionFly  wing  allows  the  wing 
to  rotate  in  addition  to  flapping.  This  chapter  presents  detailed  modeling  and 
experimental  testing  of  wing  rotation  and  lift  in  the  LionFly,  shown  schematically 
in  Figure  4.1.1.  The  goal  of  this  chapter  is  to  understand  the  flapping  and  rotation 
dynamics,  and  the  lift-producing  mechanisms  in  this  device. 

4.2  Modeling 

4.2.1  Kinematics 

The  LionFly  PAV  shown  in  Figure  4.1.1  is  a  SUEX  dry  film  compliant  slider-rocker 
mechanism  coupled  to  a  PZT-5H  bimorph  actuator.  Figure  4.2.1  shows  the  psuedo 
rigid  body  approximation  of  the  slider-rocker  mechanism.  The  flexure  hinges  are 
assumed  to  be  ideal  pins  with  torsional  stiffness  with  zero  off  axis  compliance. 
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Figure  4.1.1.  Conceptual  drawing  of  the  Penn  State  LionFly 

The  kinematics  map  the  voltage-induced  bimorph  displacement  to  flapping  angle. 
The  mechanism  is  defined  by  the  flexure  lengths,  La,  Lb,  Lq,  link  lengths,  Lab, 
Lqa,  vertical  offset,  r4,  and  the  starting  position  of  the  actuator,  xref.  The  flexure 
lengths  and  mechanism  parameters  must  satisfy 


Lab  —  Lqa  ih  +  (La  +  Lb(  1  —  7)  ~~  -^07)  >  0, 


(4.1) 
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Figure  4.2.1.  Schematic  of  slider  rocker  mechanism 


where  7  is  the  characteristic  radius  factor  of  a  long  flexure  [70].  The  kinematics 
relating  the  bimorph  displacement  to  the  flapping  angle  are 


(f>{x)  =  7 r  —  arccos 


r\  “  r3  +  rj  +  (rlo  +  x )2  \ 
2r  2  \Jr\  +  (ri0  +  xf  ) 


—  arctan 


-V). 

rlo+xJ 


(4.2) 


where  the  link  lengths  are 


77 


0 


(4.3) 


r2  =  Lqa  +  lL0  -  (1  -  7 )la, 


(4.4) 


r3  —  Lab  +  lLA  +  (1  —  ^i)LB. 


(4.5) 


The  initial  0  angle  is  determined  by  finding  the  xref  such  that  the  equilibrium 
position  xo  yields  (po{xo)  =  4>init ■  The  potential  energy  of  the  single  degree  of 
freedom  mechanism  is 


U 


k0(p(x)  T  ^  ^  kiOi^x')  T  cjkactix  37e/)  , 


i=A,B 


(4.6) 
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where  ko,  k^,  and  are  flexure  hinge  stiffnesses  of  the  three  joints  of  the  slider 
rocker  mechanism  at  0,  A,  and  B,  respectively,  kact  is  the  actuator  stiffness,  9  a 
and  9b  are  the  angles  of  the  flexure  hinges  at  A  and  B,  respectively.  The  flexure 
stiffnesses  are  calculated  by 


k 


eq 


IKq 


EfWftj 

12 Lf  ’ 


(4.7) 


where  Kq  is  a  stiffness  coefficient  dependent  on  the  range  of  motion,  Ef  is  the 
Young’s  modulus,  Wf  is  the  flexure  width,  tf  is  the  flexure  thickness,  and  Lf  is  the 
flexure  length  [70].  The  virtual  work  of  the  bimorph  is 


SW  =  GactV5x, 


(4.8) 


where  Gact  is  the  blocking  force  gain  of  the  actuator.  The  equilibrium  equation  is 


G  adV  kact(%ref  f 0? 


(4.9) 


where 


fs(x)  =  k0(p(x)^-  +  ^  kidi(x)d^ 


i=A,B 


dx 


(4.10) 


Solving  Equation  4.9  when  V  =  0  for  xref  with  x  =  xeq  =  0_1(</>o)  yields  the 
position  of  the  bimorph  needed  for  any  desired  initial  angle.  At  this  equilibrium 
position,  the  kinematics  can  be  linearized  to 


0  = 


d<j>(x) 

dx 


x  =  ax, 


X —Xeq 


(4.11) 


using  Equation  4.2. 
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Figure  4.2.2.  Schematic  of  LionFly  wing  defining  fixed  frame  xy  and  wing  frame  ry' 
and  relevant  aerodynamic  parameters 


4.2.2  Vibration 

Figure  4.2.2  shows  the  LionFly  wing  and  spar  with  two  degrees  of  freedom:  rotation 
about  the  y  axis,  0,  and  rotation  about  the  x  axis,  ip.  Using  the  kinematics 
developed  in  the  previous  section,  the  potential  energy  of  the  mechanism  is 

U  =  ^k0(p(x)2  +  ^  ki°i(x)2 

i=A,B 

^ref)  ? 


(4.12) 
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where  is  the  rotational  flexure  hinge  stiffness.  Gravitational  potential  energy 
is  negligible  compared  to  the  strain  potential  energy  stored  in  the  hinges  and 
actuator.  The  virtual  work  associated  with  the  bimorph  actuator  input  force, 
aerodynamic  forces,  aerodynamic  damping,  and  linear  viscous  damping  are 


5W  =  GactVSx  -  FN8fCp  -  -  c^Scj)  - 

=  ( Gacta~lV  -  FNrcPx  cos  0  cos  0  -  Sep 
+  (~FNrcPy  —  —  c^if)  8^ 


(4.13) 


where  rcp  is  the  distance  from  the  origin,  O,  to  the  center  of  pressure  on  the  wing, 
and  and  c ^  are  the  linear  viscous  damping  coefficients  of  0  and  0,  respectively. 
The  aerodynamic  force  normal  to  the  wing,  Fjy,  the  aerodynamic  damping  term 
associated  with  rotation,  M,/,,  and  the  center  of  pressure  location,  rcpx  and  rcpy, 
will  be  described  in  detail  in  the  next  section.  Using  the  generalized  coordinates, 
q  =  [0, 0] ,  the  kinetic  energy  is 


T  =  -qTJq 


with 


J  = 

Jyy  J s  H-  J arriyy 

J xy  J amXy 

i 

4- 

J'lfttf) 

J xy  J amXy 

J. tcx  J amxx 

J 

J'lj; 

(4,14) 


(4.15) 


where  Js  is  the  inertia  of  the  wing  spar,  and  JXXr  Jxy,  and  Jyy  are  components  of 
the  wing  inertia  in  the  wing  fixed  coordinate  system  with  associated  added  inertia 
terms  Jamxx ,  Jamxy  j  and  Jam  .  The  amplification  factor  is  large  (a  ~  3200  rad/m), 
so  the  reflected  equivalent  mass  of  the  actuator,  0.23 rnacta^2,  is  neglible  compared 
to  the  wing  and  spar  inertias.  Using  Lagrange’s  equations,  with  potential  energy 
in  Equation  4.12  linearized  about  q0  =  [axeq,  0],  kinetic  energy  in  Equation  4.14, 
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and  virtual  work  in  Equation  4.13,  the  equations  of  motion  are 


Jq  +  Cq  +  Kq  =  BQ 


where 


C 

K 

B 

Q 


c<j,  0 

1 

0  cp 

kef,  0 
0  kijj 

G actor1  —rCpx  cos(^)  cos (0)  0 

o  ~rcPy  -1 

V  Fn 


(4.16) 

(4.17) 

(4.18) 

(4.19) 

(4.20) 


where 


-  1  d2jj 
^  a2  dx2 


X=Xeq 


In  vacuum,  we  can  assume  that  the  added  mass  terms,  aerodynamic  force,  F,y, 
and  moment,  are  zero,  and  we  can  take  the  Laplace  transform  of  Eq.  4.16  to 
End  the  transfer  functions  from  applied  voltage, 


<^>(s)  _  ((G!act/'^</>)('S2  +  (cp/</p)<S  +  (/Cp/Jp))) 

_  W) 

*(*)  (g  act  J 4>ip/  s 

W)=  W) 


(4.21) 

(4.22) 


where 


D(s)  —  (s2  +  [c^j Jp)s  +  ( k J'lpj)  (s2  +  (c^/ J<p)s  +  [kfj,/ J </>)')  —  Jp) 

(4.23) 
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4.2.3  Aerodynamics 

The  aerodynamic  forces  and  moments  presented  here  use  quasi-steady  approxima¬ 
tions  and  the  blade  element  method  (BEM)  [33].  We  assume  a  hover  configuration 
with  the  device  mounted  vertically,  flapping  stroke  is  limited  to  the  horizontal 
plane,  and  zero  ambient  air  velocity.  Aerodynamic  forces  on  the  spar  are  neglected 
because  the  area  and  tip  velocity  are  significantly  smaller  than  the  wing.  Although 
many  insect  wings  deform  during  the  flap  cycle  [76],  we  assume  the  wing  spar  and 
aerofoil  are  rigid. 

4. 2. 3.1  Aerodynamic  Forces 

Figure  4.2.2  shows  relevant  aerodynamic  parameters  for  the  LionFly  wing.  The 
origin,  O,  is  defined  as  the  intersection  of  the  flapping  axis,  y,  and  the  rotation 
axis,  x.  The  wing  root  origin,  O',  is  placed  at  the  most  proximal  point  along  the 
leading  edge  of  the  aerofoil  [33],  The  axes  r  and  y'  are  shifted  from  the  x  and  y  axes 
by  xr  and  yr,  respectively.  The  position  along  the  span  of  the  wing  is  r,  and  the 
chord  along  the  span  of  the  wing  is  c(r).  The  spanwise  wing  coordinates  are  non- 
dimensionalized  by  the  span  from  wing  root  (origin  O')  to  wing  tip,  R  [77].  The 
mean  chord,  c—  A/R ,  where  A  is  the  total  wing  area,  is  used  to  nondimensionalize 
chordwise  coordinates.  The  spanwise  position  of  the  center  of  pressure  is  found  by 
finding  the  second  radius  moment  in  the  ry'  frame, 


(4.24) 


and  the  spanwise  center  of  pressure  can  be  found  by  rcpx  =  Rf2  [77] .  The  chordwise 
position  of  the  center  of  pressure  depends  on  the  instantaneous  angle  of  attack  and 
can  change  position  during  the  wing  cycle.  In  the  ry'  frame,  the  location  of  the 
center  of  pressure  for  each  blade  element  as  a  function  of  angle  of  attack,  a  is  found 
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experimentally  from  a  scaled  model  Drosophila  wing  [78], 


y'cp(r,at)  =  -c(f) 


0.82, 


led  +  0.05 


7T 


(4.25) 


In  the  xy  frame,  accounting  for  the  wing  hinge,  and  leading  edge  profile,  the 
nondimensional  chordwise  center  of  pressure  is 


Vcp 


fo  jjjr  +  VLEjf)  +  y'cp(r,  a)]  (f  +  xr)2c(r)dr 

Jq  (f  +  xr )2  c{f)dr 


(4.26) 


and  the  chordwise  center  of  pressure  can  be  found  by  rcpy  =  ycpc.  The  aerodynamic 
force  normal  to  the  wing  is 


FN  = -pair(j)\<j)\CN(a)cR3  /  (r  +  xr)2  c(r)df, 


(4.27) 


where  pair  is  the  density  of  air,  and  Cat  (ck)  is  the  aerodynamic  force  coefficient 
normal  to  the  wing,  and  a  is  the  angle  of  attack  [33] .  The  coefficients  of  lift  and 
drag  are  given  by 


CL(a)  =  1.8sin(2a),  (4.28) 

CD(a)  =  1.9  -  1.5  cos(2a),  (4.29) 

from  which  the  aerodynamic  force  coefficient  normal  to  the  wing  can  be  found  from 
the  rotational  transformation, 

Cm{cx)  =  Cl  cos(ck)  +  Co  sin(a)  3.4sin(ct),  (4.30) 

assuming  the  transverse  force  coefficient,  Ct  is  negliblc  for  a  thin  flat  plate.  In 
the  steady  case  of  hovering,  where  the  flapping  stroke  is  horizontal,  the  ‘geometric’ 
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angle  of  attack,  relative  to  the  free-stream  velocity,  can  be  used  [32]: 

a  =  ^--0-  (4.31) 

The  lift  force  is  found  using  the  rotation  angle,  ip, 

Fl  =  Fn  sin(V0  (4.32) 

4. 2. 3. 2  Aerodynamic  Moments 

The  relative  air  velocity  due  to  rotation  increases  linearly  away  from  the  hinge 
axis,  so  the  rotational  damping  is 

Mjp  =  —^Pair'lP\'lP\Crdci  Ryrd,  (4.33) 

where  the  leading  edge  position,  y0,  and  trailing  edge  position,  y1,  and  effective 
moment  arm  yrd  are 


Vo  =  yr  +  VLEpr)  -  c(r), 

(4.34) 

Vl  =  Vr  +  VLEif), 

(4.35) 

If1 

Vrd  =  1  J  [|yi|(yi)3  -  l2/o|(?/o)3]  dr 

(4.36) 

and  Crd  is  the  rotational  damping  coefficient. 


4. 2. 3. 3  Added  Mass 

Using  added-mass  coefficients  for  a  thin  fiat  wing  section[79,  33],  the  added  inertia 
terms  are 


Jyyam 


:: 


I  T  dmam  —  PairC  R  I 
' 0  4  Jo 


(f  +  xr)2  c(r)2  dr , 


(4.37) 
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Figure  4.3.1.  Photograph  of  fabricated  prototype  LionFly  LF1203R20 


-R 


JxUam 


7T 


ryhdmam  =  -pairc2R2  /  (r  +  xr)yhc(r)2df, 


rR 


7T 


Jxxam  =  (yh  +  7^c(r)  )dmam  =  -pairc  R  /  j/fcc(r)  dr, 


(4.38) 

(4.39) 


where  y h  is  the  nondimensional  distance  from  the  midpoint  to  the  rotational  axis, 


Vh{r)  =  g c{f )  -  yLE(r)  ~  yr- 


(4.40) 


4.3  Experimental  Testing 


Figure  4.3.1  shows  the  LionFly  mechanism  fabricated  using  the  process  developed 
in  Section  3.2.2  and  mounted  onto  a  glass  slide  for  characterization.  A  2  mm  x 
16  mm  PZT-5H  bimorph  actuator  is  mounted  using  precision  linear  stages.  Using 
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Figure  4.3.2.  Photograph  of  experimental  setup  used  for  high  voltage  experiments 


the  stages,  the  wing  is  positioned  at  an  initial  angle  0o  and  fixed  by  clamping  the 
bimorph  with  a  2  mm  x  2  mm  glass  square  bonded  with  cyanoacrylate.  Gold 
reflective  surfaces  are  attached  to  certain  positions  on  the  wing  for  small  signal  laser 
vibrometer  measurement.  A  Polytec®  laser  vibrometer  measures  the  displacement 
at  the  test  points  on  the  wing  surface..  The  vibrometer  uses  a  10 x  objective  lens 
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Figure  4.3.3.  Photograph  of  experimental  setup  used  for  lift  measurement  experiments 


to  focus  the  laser  to  a  25  /im  point  and  returns  displacement  and  velocity  data  to 
LabView®  with  sub  micron  precision.  For  small  signal  actuation,  the  LabView® 
DAQ  system  can  provide  a  5  V  amplitude  signal,  but  for  high  voltage  experiments 
(25  —  75  V  amplitude  with  75  V  DC  bias),  a  Matsusada®  AP-1B3  amplifies  the 
DAQ  output.  Figure  4.3.2  shows  a  photograph  of  the  experimental  setup  used  to 
measure  high  voltage  wing  trajectories  using  stroboscopic  photography.  LabView 
controls  the  digital  SLR  camera,  strobe  light,  and  the  Matsusada®  amplifier.  A 
frequency  sweep  is  conducted  by  adjusting  the  strobe  frequency  with  the  actuation 
frequency  and  taking  twenty  pictures  throughout  the  wing  flapping  cycle  at  each 
frequency  point.  The  photographs  are  processed  in  custom  Matlab®  software  and 
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wing  (p{t)  and  ip(t)  data  are  extracted. 

Figure  4.3.3  shows  a  photograph  of  the  lift  measurement  setup.  The  LionFly 
is  positioned  vertically  and  bolted  to  an  aluminum  lever  arm  with  a  razor  blade 
fulcrum  that  rests  on  an  Aurora  Scientific®  Inc.  402A  force  transducer  with  10 
yN  resolution  and  2kHz  bandwidth.  The  lever  constrains  the  force  measurement 
to  a  single  degree  of  freedom,  preloads  the  transducer,  and  allows  dynamic  loads  to 
be  transmitted  to  the  transducer.  When  the  device  is  actuated,  the  lift  generated 
decreases  the  measured  load  on  the  transducer  proportional  to  the  leverage,  and 
average  lift  data  is  taken  from  100  cycles. 

Flexure  and  membrane  thicknesses  are  measured  using  a  Zygo®  NewView  100C 
optical  profilometer  and  vein  thickness  using  a  micrometer.  The  initial  angle  and 
vertical  offset  were  measured  using  high  resolution  macro  photographs  with  cali¬ 
brated  pixels.  A  three  dimensional  model  of  the  wing,  spar,  veins,  and  test  points 
was  created  in  Matlab  using  discretized  images  and  the  measured  vein  and  mem¬ 
brane  thickness,  tvein  and  tmem,  respectively.  The  inertia  terms  were  calculated  by 
numerically  integrating  the  model  wing  using  the  reported  value  for  SUEX  density, 
Psuex •  As  a  check,  the  mass  was  experimentally  measured  and  verified  to  match 
the  calculated  mass  of  the  wing.  Using  the  measured  value  of  wing  span,  R ,  and 
the  profile  of  the  wing  and  spar,  all  aerodynamic  parameters  can  be  numerically 
calculated. 

4.3.1  Experimental  Results 

Out-of-plane  laser  vibrometer  data  is  transformed  into  flapping  and  rotation  angles 
using  the  displacement  of  test  points  on  the  wing,  shown  in  Figure  4.3.1 


rtPi  =  Xii  +  ytj  +  Zik, 


(4.41) 
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Figure  4.3.4.  Flapping  frequency  response  4>(s)  of  LionFly  LF1203  (0  VDC/5  VAC): 
Experimental  data  in  air  (blue  circles),  theoretical  response  in  air  (blue  solid),  experi¬ 
mental  data  in  vacuum  (green  diamonds),  theoretical  response  in  vacuum  (green  dashed), 
simulated  response  in  air  (red  dashed  dot),  simulated  response  in  vacuum  (black  dotted). 


where  Xi  and  r/i  are  the  x  and  y  distance  from  the  0  and  0  axes,  respectively,  and  zt 
is  the  measured  vibrometer  displacement,  given  by  Zi  =  5mease:’ev^9s ,  where  Smeas 
is  a  vector  of  magnitudes  of  the  displacement  frequency  response,  6y  is  the  phase 
of  the  voltage  signal,  and  9s  is  the  phase  of  the  displacement  signal.  The  vector 
perpendicular  to  the  wing  plane  can  be  written  as 

n  =  (rtp 3  -  rtp4 )  x  (rtps  -  ftpA)  =  xni  +  ynj  +  znk. 


(4.42) 
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Figure  4.3.5.  Rotation  frequency  response  'k(s)  of  LionFly  LF1203  (0  VDC/5  VAC): 
Experimental  data  in  air  (blue  circles),  theoretical  response  in  air  (blue  solid),  experi¬ 
mental  data  in  vacuum  (green  diamonds),  theoretical  response  in  vacuum  (green  dashed), 
simulated  response  in  air  (red  dashed  dot),  simulated  response  in  vacuum  (black  dotted). 


and  the  transformed  angles  are  calculated  as 

0  =  arctan  ^  ,  (4.43) 

0  =  arctan  ^  •  (4.44) 

Figures  4.3.4  and  4.3.5  show  the  frequency  response  of  flapping  and  rotation, 
respectively,  under  applied  voltage  of  5  V  amplitude,  using  Equations  4.21  and 
4.22,  respectively,  and  also  using  the  nonlinear  aerodynamics  in  Equation  4.16. 
Tables  4.3.1-  4.3.4  show  the  selected  parameter  and  values  used  in  the  modeling. 
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Table  4.3.1.  Constant  parameters  and  values 


Constants  {Units) 

Description 

Value 

Ef  ( GPa ) 

Young’s  modulus  of  SUEX 

2.80 

7 

characteristic  radius  factor 

0.8517 

I<e 

stiffness  coefficient 

2.6762 

Psuex  (kg/m3) 

density  of  air 

1190.00 

Pair  (kg/m3) 

density  of  air 

1.28 

cN 

aerodynamic  normal  force  coefficient 

3.40 

Crd 

aerodynamic  rotational  damping  coefficient 

5.00 

Table  4.3.2.  Selected  LionFly  actuator  parameters  and  values 


Parameter 

Description 

Value 

Gact. 

blocking  force  gain 

328.29 

pN/V 

act 

stiffness 

658.83 

N/m 

act 

mass 

81.69 

mg 

Xref 

reference  position 

53.40 

pm 

The  linear  viscous  damping  terms  are  tuned  to  match  the  experiment  in  air  and  in 
vacuum.  All  other  parameters  are  independently  measured  or  calculated  as  stated 
earlier.  The  results  show  excellent  agreement  for  the  first  mode,  near  38  Hz.  The 
response  in  vacuum  is  much  larger  than  in  air  due  to  aerodynamic  damping.  The 
first  flap  and  rotation  natural  frequencies  are  also  reduced  from  vacuum  to  air, 
showing  the  effect  of  added  aerodynamic  inertia.  The  simulated  and  experimen¬ 
tal  results  have  the  same  decrease  in  natural  frequency  so  the  model  accurately 
captures  the  added  mass  effect.  The  flap  and  rotation  amplitudes  and  phases  are 
roughly  equal  in  the  first  mode,  indicating  strong  coupling  and  potentially  poor 
lift  production,  respectively.  At  this  low  level  of  excitation  the  aerodynamic  forces 
are  not  sufficient  to  rotate  the  wing  and  change  the  zero  degree  relative  phase 
between  flap  and  rotation.  It  is  therefore  expected  that  the  lift  produced  by  this 
wing  motion  will  be  extremely  low  because  the  phase  difference  is  far  from  the 
optimal  value  of  90  degrees  [11].  The  frequency  response  for  the  nonlinear  sim- 
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Table  4.3.3.  LionFly  mechanism  parameters  and  values 


Parameter 

Description 

Value 

Lab 

length  of  coupler  link 

3.25 

(mm) 

Loa 

length  of  rocker  link 

0.75 

(mm) 

La 

length  of  flexure  hinge  A 

1.00 

(mm) 

Lb 

length  of  flexure  hinge  B 

0.25 

(mm) 

Lo 

length  of  flexure  hinge  0 

1.00 

(mm) 

r. 4 

vertical  offset  of  slider 

0.77 

(mm) 

tf 

mean  flexure  thickness 

24.7 

(/im) 

kA 

stiffness  of  flexure  hinge  A 

18.43 

(fiN  ■  m) 

bn 

stiffness  of  flexure  hinge  B 

73.73 

(fj,N  ■  m) 

k0 

stiffness  of  flexure  hinge  0 

44.88 

(fiN  ■  m) 

4>o 

initial  0  angle 

2.50 

(deg.) 

Wlmech 

mass 

21.62 

(mg) 

ulation  also  matches  the  experimental  results  at  low  voltage.  At  each  frequency 
point,  the  nonlinear  simulation  is  run  to  steady  state  and  the  amplitude  and  phase 
calculated  numerically.  The  vacuum  results  show  that  the  nonlinear,  linear,  and 
experimental  results  all  agree  for  the  first  mode,  matching  both  magnitude  and 
phase  for  both  flap  and  rotation.  The  simulated  results  in  air  do  not  include  the 
tuned  aerodynamic  damping  of  the  linear  model  but  match  the  linear  and  exper¬ 
imental  data  very  well  nonetheless.  This  validates  the  accuracy  of  the  nonlinear 
aerodynamic  model  at  low  amplitudes. 

The  agreement  between  models  and  experiment  at  the  second  mode,  near 
135  Hz  is  also  reasonable.  The  flap  and  rotation  amplitudes  and  frequencies  are 
slightly  larger  in  vacuum  than  in  air,  again  indicating  the  effects  of  aerodynamic 
drag  and  added  mass.  The  model  slightly  over  predicts  the  second  natural  fre¬ 
quency,  probably  due  to  unmodeled  compliance  in  the  LF1203  wing  and  mecha¬ 
nism.  The  second  mode  couples  flap  and  rotation  but  rotation  has  a  substantially 
larger  amplitude  relative  to  the  first  mode.  Thus,  the  first  mode  is  primarily  a  flap 
mode  and  the  second  is  primarily  a  rotation  mode  but  both  are  strongly  coupled 
and  cannot  be  accurately  analyzed  or  understood  with  a  decoupled  model.  Decou- 
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Table  4.3.4.  LionFly  wing  parameters  and  values 


Parameter 

Description 

Value 

R 

wing  span  from  O'  to  tip 

20.10 

(mm) 

c 

mean  chord 

5.82 

(mm) 

A 

area  of  wing 

116.89 

(mm2) 

xr 

x  distance  from  0  to  O' 

3.00 

(mm) 

Ur 

y  distance  from  0  to  O' 

-0.21 

(mm) 

rcPx 

x  distance  from  0  to  CP 

13.43 

(mm) 

rtp3 

(x,  y )  distance  from  0  to  ftp3 

(10.47,-1.88) 

(mm) 

rtP4 

(x,  y)  distance  from  0  to  ftp4 

(11.60,-7.86) 

(mm) 

nP5 

(x,  y)  distance  from  0  to  ftp5 

(22.87,-3.29) 

(mm) 

ttp 

effective  thickness  of  PET  test  surface 

75.04 

(ym) 

tvein 

average  thickness  of  veins 

250 

(ym) 

tmem 

average  thickness  of  membrane 

250 

(ym) 

tip 

flexure  thickness  of  rotational  hinge 

25.74 

(ym) 

stiffness  of  wing  rotational  hinge 

44.05 

(yN  ■  m) 

mw 

mass  of  wing 

10.44 

(mg) 

ms 

mass  of  spar 

8.62 

(mg) 

Js 

spar  inertia  about  y  axis 

170.82 

(g  ■  mm2) 

Js 

° am 

added  spar  inertia  about  y  axis 

1.40 

(g  ■  mm2) 

Jyy 

wing  inertia  about  y  axis 

1738.17 

(g  ■  mm2) 

'kyya.m 

added  inertia  about  y  axis 

117.84 

(g  ■  mm2) 

J xy 

product  of  inertia  about  y  axis 

-461.07 

(g  ■  mm2) 

Jxyam 

product  of  added  inertia  about  y  axis 

-37.29 

(g  ■  mm2) 

J XX 

wing  inertia  about  x  axis 

175.17 

(g  ■  mm2) 

JxXam 

added  inertia  about  x  axis 

14.14 

(g  ■  mm2) 

C(j),vac 

linear  viscous  damping  of  (p  in  vac 

10.02 

(mN  ■  ym  ■  s) 

C(f),air 

linear  viscous  damping  of  (p  in  air 

19.60 

(mN  ■  ym  ■  s) 

C"i{),vac 

linear  viscous  damping  of  ip  in  vac 

4.46 

(mN  ■  ym  ■  s) 

C"il),air 

linear  viscous  damping  of  ip  in  air 

5.70 

(mN  ■  ym  ■  s) 

pled  analyses  may  be  appropriate  for  systems  with  light  wings  where  the  rotation 
resonance  is  much  higher  than  the  flap  resonance  (10  times)  but  the  stiffness  is  low 
enough  to  allow  rotation  under  aerodynamic  loading  at  the  flap  resonance.  This 
is  clearly  not  the  case  for  LF1203.  Although  the  two  resonances  are  separated  by 
a  factor  of  three,  both  modes  involve  significant  flap  and  rotation. 
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Figure  4.3.6.  Stroboscopically  photographed,  image  processed,  and  fitted  flapping,  </>(t) 
and  rotation  angles,  ip(t),  and  calculated  lift  force  (75  VDC / 75  VAC)  and  36  Hz:  (a) 
Experimental  4>(t)  (circles)  and  ip(t)  (diamonds)  and  fitted  4>(t)  (solid)  and  ip(t)  (dashed), 
and  (b)  calculated  lift  force  average  (green  dashed)  and  calculated  using  experimental 
( f>(t )  and  -0(t)  data  (blue  solid). 


At  higher  actuation  voltages,  the  laser  vibrometer  cannot  be  used  for  angle 
measurement,  so  stroboscopic  photography  and  image  processing  measure  the  flap 
and  rotation  angles  versus  time.  The  wing  trajectory  points  measured  from  pho¬ 
tographs  are  linearly  interpolated  and  fitted  using  the  magnitude  and  phase  from 
the  highest  four  frequency  components  as  calculated  using  the  fast  fonrier  trans¬ 
form.  Figure  4.3.6(a)  shows  a  time  trace  of  one  period  of  measured  and  fitted  wing 
trajectories  at  36  Hz.  Figure  4.3.6(b)  shows  the  lift  calculated  from  the  measured 
angle  trajectories  using  Equations  4.27,  4.30,  and  4.32.  While  the  lift  varies  widely 
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Figure  4.3.7.  Experimental  and  simulated  flapping,  rotation,  and  rela¬ 

tive  phase,  Z< p(t )  —  Aip(t),  versus  frequency:  (a)  experimental  (p(t)  (circles)  and  ip(t) 
(diamonds)  amplitudes  and  simulated  4>(t)  (solid)  and  ip(t)  (dashed)  amplitudes  at  75 
VDC/ 75  VAC  (blue),  and  75  VDC/ 25  VAC  (green),  and  (b)  experimental  relative 
phase  (circles)  and  simulated  relative  phase  (solid)  at  75  VDC /75  VAC  (blue)  and  75 
VDC/25  VAC  (green).  Data  points  corresponding  to  Figure  4.3.6(a)  are  filled. 


during  the  cycle,  the  mean  value  is  a  small  11.52  fiN.  The  reason  for  the  low  lift 
is  that  the  </>(t)  and  time  trajectories  are  almost  completely  in  phase  at  this 
frequency.  The  trajectory  lags  the  (j){t)  trajectory  by  approximately  12.06 
deg.,  much  less  than  the  optimal  value  of  90  degrees  that  produces  maximal  lift. 

Figure  4.3.7  shows  the  measured  and  simulated  high  voltage  wing  trajectories 
and  relative  phase  difference  between  flapping  and  rotation  under  applied  voltage  of 
75  V  DC  bias,  and  amplitudes  of  25  V  (75  VDC/25  VAC )  and  75  V  (75  V DC /75 
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VAC)  versus  frequency.  The  23  deg.  amplitude  </>(£),  22  deg.  amplitude  ^(t),  and 
12  deg.  relative  phase  results  from  Figure  4.3.6(a)  are  shown  as  filled  symbols  at 
37  Hz.  At  25  V,  the  peak  wing  angles  occur  at  a  resonance  slightly  higher  than 
the  low  excitation  resonance  because  the  75  V  offset  shifts  the  linearization  point, 
slightly  decreasing  amplification,  and  thereby  increasing  resonance  frequency.  The 
nonlinearity  of  the  mechanism  is  apparent  in  the  shift  in  resonance  frequency  when 
amplitude  is  again  increased  in  the  (75V DC /75V AC)  case.  In  both  cases,  flapping 
and  rotation  magnitudes  are  well  predicted,  however  the  phase  differences  are 
under  predicted.  The  phase  difference  is  highly  dependent  on  the  aerodynamic 
forces  acting  upon  the  wing,  in  which  changes  in  center  of  pressure  location  and 
uncertainties  in  aerodynamic  force  coefficients  contribute  to  error. 

Figure  4.3.8  shows  the  average  lift  frequency  response  of  the  LionFly  using 
experimentally  measured  lift,  calculated  lift  from  measured  wing  trajectories,  and 
simulated  lift.  The  11.52  gN  average  lift  result  from  Figure  4.3.6(b)  is  shown  as  the 
filled  symbol  at  36  Hz.  The  experimentally  measured  lift  contains  data  from  three 
runs  at  each  applied  voltage,  showing  consistent  performance  of  the  lift  measure¬ 
ment  system.  The  measured  lift  is  larger  than  predicted  by  the  calculated  lift  from 
wing  trajectory  data  and  simulation.  Accuracy  of  measured  lift  depends  strongly 
on  the  relative  phase  difference  in  the  measured  wing  trajectory  data,  which  is 
limited  by  temporal  and  pixel  resolution.  The  Reynold’s  number  approximated 
by  Re  =  (2 <f>Rc)/is,  where  v  is  the  kinematic  viscosity  of  air,  v  «  20  —  26 mm2/s, 
yields  ~  1500,  higher  than  the  Reynold’s  number  («  100)  used  in  finding  aerody¬ 
namic  coefficients  in  [11].  Additionally,  unmodeled  aerodynamic  phenomena,  such 
as  rotational  lift,  wake  capture,  etc.,  could  be  responsible  for  these  differences. 
Simulated  lift  values  suffer  from  low  magnitude  as  expected  from  the  underpre¬ 
dicted  phase  data.  Both  the  25  V  case  and  the  75  V  case  show  maximum  lift  at  a 
frequency  slightly  higher  than  the  resonance  observed  in  the  wing  trajectory  data, 
due  to  increased  relative  phase  difference  with  slightly  off-resonance  flapping  and 
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Figure  4.3.8.  Average  value  lift  versus  frequency:  experimentally  measured  lift  (blue), 
calculated  lift  from  wing  measurement  (green)  at  75  VDC  / 75  VAC  (circles),  75  VDC / 25 
VAC  (diamonds),  and  simulated  lift  (red)  at  75  VDC /75  VAC  (solid)  and  75  VDC  1 25 
VAC  (dashed).  Data  point  corresponding  to  Figure  4.3.6  is  filled. 


rotation  amplitudes.  Despite  uncertainties  in  aerodynamic  parameters,  the  mea¬ 
sured  wing  trajectories  and  lift  validate  the  key  criteria  required  for  lift  production: 
large  flapping,  rotation,  and  relative  phase  trajectories. 


4.4  Model-Based  Redesign 

The  lift  generated  by  LF1203  is  much  smaller  than  what  would  be  required  to  hover 
(1195  /J.N),  not  including  the  payload  weight.  The  vehicle  design  is  constrained  by 
the  current  fabrication  process  that  produces  a  minimum  SUEX  thickness  of  20  /im. 
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Figure  4.4.1.  Simulated  wing  trajectories  and  lift  force  versus  time  of  improved  device 
of  the  redesigned  LionFly  at  applied  voltage  of  (0  VDC/  300  VAC)  and  123  Hz:  (a) 
flapping,  4>(t)  (blue  solid),  and  rotation,  ijj(t)  (green  dashed),  (b)  lift  force,  Fjj  (blue 
solid),  and  average  lift  (green  dashed). 


This  minimum  thickness  is  important  because  it  governs  the  stiffness  of  the  rotation 
hinge  and  the  mass  of  the  wing.  Optimal  lift  production  requires  large  wing  flap 
and  rotation  and  a  phase  difference  of  around  90  degrees.  A  thinner  rotation  hinge 
would  allow  larger  rotation  but  lower  the  “rotation  natural  frequency,  amplifying 
the  modal  coupling.  The  modal  solutions  typically  do  not  provide  good  phase 
difference  so  this  coupling  must  be  overcome  by  the  aerodynamics.  To  maintain 
a  high  rotation  frequency,  the  mass  of  the  wing  must  also  be  reduced.  The  mass 
of  the  wing,  however,  also  dominates  the  flap  natural  frequency  so  it  is  not  clear 
that  the  spectral  separation  between  the  two  modes  would  be  improved.  In  this 
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section,  the  LionFly  is  redesigned  using  the  model  to  produce  lift  equal  to  its 
weight,  enabling  a  hovering  device  in  theory.  This  is  by  no  means  an  optimized 
device  and  the  previous  section  showed  underprediction  of  lift  by  the  model  so 
the  actual  performance  in  practice  may  be  better  than  predicted.  The  purpose  of 
this  section,  however,  is  to  demonstrate  that  if  the  constraints  on  the  fabrication 
process  are  removed  then  a  higher  performing  PAV  becomes  a  possibility. 

From  the  experiment,  it  is  clear  that  the  wing  rotation  of  LF1203  is  nearly  in 
phase  with  flapping  and  aerodynamic  forces  are  small  compared  to  the  wing  inertia. 
If  the  wing  inertia  is  reduced,  resonance  frequencies  are  increased,  potentially 
increasing  lift  production.  Fabrication  of  a  wing  with  vein  thickness  of  30  gm 
and  membrane  thickness  of  1  gm  (much  smaller  than  the  250  gm  and  20  gm 
contraints  of  the  LF1203  fabrication  process),  wing  mass  can  be  reduced  from 
10.4  mg  to  1.35  mg,  reducing  the  inertia  accordingly.  A  second  design  change  is 
the  specification  of  an  initial  flap  angle  of  45  deg.  This  is  difficult  to  dial  in  exactly 
using  the  current  process  so  DC  voltages  are  applied  to  adjust  the  initial  flap  angle. 
Specifying  an  initial  flap  angle  of  45  deg  provides  a  wide  operating  range  and 
lower  amplification  factor.  The  reduced  amplification  factor  raises  the  resonance 
frequencies,  a  positive  result,  but  also  increases  the  required  voltage  to  obtain  the 
desired  flap  amplitude,  a  negative  result.  The  overall  voltage  amplitude  is  therefore 
increased  to  300V,  within  the  operating  range  of  the  LF1203  bimorph.  The  final 
design  change  is  a  reduction  in  the  rotational  flexure  thickness,  £</,  =  7.5  fim,  and 
the  flexure  length,  =  352  gm.  These  three  design  changes  greatly  improve  the 
amplitudes  and  relative  phase,  potentially  leading  to  higher  lift  production. 

Figure  4.4.1  shows  the  simulated  flapping  and  rotation  wing  trajectories  versus 
time  at  the  maximum  average  lift  frequency  of  the  redesigned  LionFly  (123  Hz). 
Compared  with  Figure  4.3.6,  significant  phase  shift  between  flapping  and  rotation 
combined  with  large  rotation  angles  should  lead  to  much  larger  lift  forces.  The 
increased  resonance  frequency  also  quadratically  increases  aerodynamic  forces. 
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Figure  4.4.2.  Simulated  flapping,  0(f),  rotation,  0(f),  and  relative  phase,  Z0(t)— Z0(f), 
versus  frequency  for  the  redesigned  LionFly  (0  VDC /  300  VAC  -  blue  and  0  VDC /  150 
VAC  -  green):  (a)  0(f)  (solid)  and  0(f)  (dashed),  and  (b)  relative  phase.  The  flapping 
resonance  frequency  is  highlighted  at  117  Hz  (gray  shaded). 


Figure  4.4.2  shows  the  simulated  high  voltage  frequency  reponse  of  the  re¬ 
designed  LionFly  using  applied  voltages  of  150  V AC  and  300  VAC.  These  results 
show  much  larger  rotation,  similar  flapping,  and  most  importantly,  the  larger  rel¬ 
ative  phase  compared  to  LF1203. 

Figure  4.4.3  shows  the  simulated  lift  frequency  response,  with  peak  lift  at  1250 
/j,N.  The  peak  lift  occurs  at  a  frequency  slightly  above  the  flapping  resonance, 
similar  to  the  experiment.  This  effect  increases  as  a  function  of  voltage  ampli¬ 
tude.  The  simulated  lift  brings  the  LionFly  within  a  lift  to  weight  ratio  of  one, 
theoretically  allowing  the  device  to  hover. 
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Figure  4.4.3.  Simulated  average  steady  state  lift  versus  frequency  of  the  redesigned  Li- 
onFly  (0  VDC /  300  VAC  -  blue  and  0  VDC /  150  VAC  -green).  The  flapping  resonance 
frequency  is  highlighted  at  117  Hz  (gray  shaded) 


Chapter  5 

Conclusions  and  Future  Work 


In  this  work,  the  design,  fabrication,  modeling,  and  experimental  validation  of  an 
NAV  scale  clapping  wing  mechanism  is  presented.  A  flexure  hinge  allows  passive 
wing  rotation  for  the  clapping  wing  mechanism  that  has  the  potential  to  increase 
thrust.  Analytical  models  of  wing  flapping  and  rotation  are  derived  and  validated 
using  experimental  wing  trajectory  results.  The  Penn  State  NAV  prototype  is 
experimentally  shown  to  provide  approximately  54o  peak  to  peak  wing  rotation, 
14o  peak  to  peak  flapping  angle,  and  0.21  mN  of  thrust  at  9.5  Hz.  At  25.5  Hz , 
the  prototype  produces  a  maximum  of  1.34  mN  of  thrust.  The  PSNAV  model 
accurately  predicts  the  wing  resonances  in  the  experimental  prototype.  Model- 
predicted  thrust  is  lower  than  the  experimentally  measured  values,  however. 

This  work  also  presents  a  simple,  low  cost,  precise,  and  rapid  process  for  mono¬ 
lithic  fabrication  of  SUEX  M  dry  film  to  prototype  PAV  mechanisms  and  wings. 
Several  LionFly  prototypes  were  fabricated  and  experimentally  tested.  Theoret¬ 
ical  and  experimental  results  have  excellent  agreement  validating  the  compliant 
mechanism  kinematics  and  aerodynamic  added  mass  and  damping.  High  voltage 
tests  show  a  flapping  angle  of  55°  at  150  V  amplitude  peak  to  peak  with  150  V  DC 
offset  at  51  Hz  resonance.  The  results  demonstrate  the  potential  of  this  process  to 
fabricate  PAV  scale  flapping  mechanisms  and  wings,  and  the  ability  of  the  rnech- 
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anism  to  amplify  the  small  displacement  of  the  piezoelectric  bimorph  into  large 
flapping  angle.  Consistent  performance  from  multiple  prototypes  demonstrate  the 
reliable  and  repeatable  nature  of  the  fabrication  process. 

Finally,  this  work  presents  detailed  modeling  and  experimental  testing  of  the 
flapping  and  rotation  dynamics,  and  the  lift  producing  mechanisms  in  the  LionFly. 
The  linear  vibration  model  augmented  with  nonlinear  aerodynamic  forces  was  val¬ 
idated  using  experimental  testing  with  a  laser  vibrometer  and  accurately  predicts 
small  amplitude  wing  dynamics  in  air  and  vacuum.  At  higher  amplitudes,  the 
model  can  sufficiently  predict  wing  trajectory  amplitudes,  but  phase  measurement 
and  simulation  have  slight  error.  The  LionFly  produces  46  deg.  flap  and  44  deg. 
rotation  amplitude  peak  to  peak  with  relative  phase  of  12  deg.,  and  maximum  lift 
of  71  gN  at  37  Hz.  Calculated  lift  values  using  measured  wing  trajectories  and 
simulated  lift  values  underpredict  measured  lift.  Sensitivity  to  phase  measurement, 
unmodeled  aerodynamic  mechanisms,  and  uncertainty  in  aerodynamic  parameters 
each  contribute  to  error.  Despite  this,  the  measured  wing  trajectories  and  lift 
validate  the  key  criteria  required  for  lift  production:  large  flapping,  rotation,  and 
relative  phase  trajectories.  By  reducing  the  inertia  of  the  wing  and  tuning  the 
rotational  hinge  stiffness,  a  redesigned  device  is  simulated  to  produce  lift  to  weight 
ratio  of  one. 


5.1  Future  Work 

5.1.1  Fabrication 

The  current  fabrication  process  has  many  advantages  in  simplicity,  time,  and  cost. 
Using  the  validated  model,  we  have  shown  that  reducing  the  membrane  thickness 
to  1  gm  along  with  wing  venation  to  30  gm  will  cause  a  significant  increase  in 
performance.  However,  the  thickness  range  for  the  thin  region  is  between  20  -  75 
gm.  The  minimum  thickness,  found  through  a  reduction  of  yield  due  to  membrane 
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a)  AZ  4562  on  Si,  blanket  exposure 
Cross  section  of  AA'  (Prefab) 


f 


A 


b)  SUEX  laminated  onto  substrate 


d)  2nd  3 10  nm  UV  dose  defines  hinges 


e)  365  nm  UV  dose  defines  links 


f)  PEB  crosslinks  exposed  SUEX 


g)  Parylene  C  deposition 


h)  1 827  Patterned  for  02  Plasma 


i)  02  Plasma  defines  membrane 


|  AZ  4562 
|  SC- 1827 
■  Si 


□  Unexposed  SUEX 
I  Exposed  SUEX 


Crosslinked  SUEX 


Parylene  C 


i)  PGMEA  develops  SUEX,  releases 
Cross  section  of  AA'  (Postfab) 


Mechanism 


Figure  5.1.1.  Potential  fabrication  process  for  PAV  flapping  wing  mechanisms:  a) 
prepare  of  Si  wafer  with  an  AZ  4562  sacrificial  layer,  b)  laminate  SUEX  onto  substrate, 
c)  define  venation  through  310  nm  UV  dose,  d)  define  flexure  hinges  through  2 nd  310 
nm  UV  dose,  e)  define  links  through  365  nm  UV  dose,  f)  crosslink  SUEX  through  post 
exposure  bake  at  65°C  for  30  minutes,  g)  1  fim  parylene  C  deposition,  h)  pattern  SC-1827 
to  mask  membrane  (spin  SC-1827,  soft  bake  at  95°C  for  2  minutes,  expose  using  310  nm 
UV  dose,  develop  in  acetone,  N2  dry),  i)  oxygen  plasma  removes  undesired  parylene  C, 
j)  develop  and  release  in  PGMEA,  dry  using  N2  and  a  50°C  vacuum  oven. 


perforation,  is  feature  size  dependent  and  is  strongly  correlated  to  stress  that 
develops  during  the  post  exposure  bake.  More  characterization  and  tests  need  to 
be  conducted  to  better  understand  this  issue.  Alternatively,  the  process  could  be 
modified  to  include  a  separate  material  for  the  1  /im  membrane  layer. 

Figure  5.1.1  shows  a  potential  fabrication  process  to  implement  the  design 
changes  proposed  in  section  4.4.  The  previous  process  only  used  photolithogra¬ 
phy  which  allowed  SUEX  to  be  used  as  a  substrate.  In  this  process,  SUEX  must 
be  laminated  to  a  substrate  for  additional  processing.  The  process  begins  with 
preparing  a  sacrificial  layer  for  releasing  the  device,  shown  in  Figure  5.1.1  a).  AZ 
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4562  has  recently  been  used  as  a  sacrificial  layer  in  SU-8  processing  [80]  and  is 
much  simpler  to  use  than  other  sacrificial  layers  such  as  aluminum  which  must  be 
deposited  by  sputtering  or  evaporation.  The  SUEX  must  be  hot  roll  laminated 
onto  the  AZ  4562/Si  substrate  shown  in  Figure  5.1.1  b).  The  vein  layer  and  the 
flexure  hinge  layer  can  both  be  defined  using  different  310  nm  UV  doses  as  shown 
in  Figure  5.1.1  c)  and  d),  respectively.  Figure  5.1.1  e)  shows  the  rigid  link  defini¬ 
tion  using  365  nm  UV  dose.  The  post  exposure  bake  in  Figure  5.1.1  f)  at  65°C  for 
30  minutes  crosslinks  the  SUEX.  Parylene  C  is  a  suitable  candidate  for  membrane 
material  because  it  is  easily  deposited  at  room  temperature  at  thicknesses  of  1 
fim  and  can  withstand  the  chemical  development  used  for  SUEX.  Figure  5.1.1  g) 
shows  the  parylene  C  deposition,  and  Figure  5.1.1  h)  shows  patterned  SC-1827. 
This  step  involves  spinning  SC-1827,  soft  baking  The  membrane  must  be  defined 
using  an  oxygen  plasma  etch,  which  requires  a  masking  layer  of  SC-1827.  The  1 
/im  parylene  layer  will  etch  rather  quickly  in  comparison  to  the  approximately  3 
/jm  layer  of  photoresist,  forming  the  wing  membrane  in  Figure  5.1.1  i).  The  final 
development  in  PGMEA  dissolves  the  uncrosslinked  SUEX  and  AZ4562,  simulta¬ 
neously  releasing  the  device.  The  device  is  dried  using  N2  and  a  vacuum  oven  at 
50° C  vacuum  oven. 

This  potential  process  assumes  SUEX  will  not  be  negatively  affected  by  the  SC- 
1827  soft  bake  temperatures,  exposure,  or  02  plasma.  The  SC-1827  data  sheet  soft 
bake  temperature  of  110°C  can  be  adjusted  to  95°C  with  extra  time  to  be  within 
reported  values  for  PEB  of  SUEX.  The  exposure  dose  required  for  SC-1827  is  an 
order  of  magnitude  smaller  than  SUEX,  and  no  post  exposure  bake  is  conducted 
subsequently  so  no  undesired  crosslinking  of  SUEX  should  occur.  Defining  the 
parylene  C  layer  using  the  02  plasma  could  require  tuning  of  SC-1827  thickness 
and  etch  time.  This  process  presents  one  possible  solution  to  utilize  the  current 
benefits  without  significant  increase  in  processing  steps. 
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Figure  5.1.2.  Parameterization  schematic  of  simplified  Lionfly  wing  design 


5.1.2  Design 


Optimizing  the  wing  under  the  current  fabrication  constraints  would  be  invaluable 
for  future  prototypes.  To  explore  the  design  space,  simplifying  the  wing,  spar,  and 
rotational  hinge  allows  parameterization  based  on  a  few  variables.  Based  on  the 
current  prototype,  a  simplified,  rectangular  design  can  be  formed  as  a  function  of 
wing  span,  i?,  and  average  chord,  c,  by  Figure  5.1.2.  For  a  constant  aspect  ratio, 
AR ,  the  average  chord  can  be  written  as,  c  =  R/AR.  The  wing  planform  has 
thickness,  tavg,  and  the  inertia  matrix  is 


J 


_  +  JD  4 

^avg1 1 


157.31  34.87 
34.87  10.76 


(5.1) 


(5.2) 


The  aerodynamic  moments  on  acting  on  0  and  0  rewritten  in  terms  of  R  are 

M 0  =  R5  ^O.18760|0|  cos2(0)  cos(0)  j  ,  (5.3) 

Mv,  =  R5  ^0.0396 cos(0)0|0|  +  O.O1170|0|)  .  (5.4) 

This  shows  that  the  inertia  terms  have  a  R4  dependence  while  the  aerodynamics 
have  a  R 5  dependence.  The  rotational  hinge  stiffness,  k^,  determined  strongly  by 
the  flexure  thickness,  t,/, ,  is  another  design  variable  which  can  be  used  in  optimiza¬ 
tion.  Selecting  an  appropriate  tavg,  ty  between  20  to  75  /im,  the  current  model  can 
be  used  to  find  an  optimal  wing  span  for  a  fixed  actuator  and  mechanism. 

Further  optimization  would  include  the  design  space  of  the  actuator  and  mech¬ 
anism  which  govern  the  flapping  degree  of  freedom  stiffness,  amplification,  a, 
bimorph  force  gain,  Gact,  and  the  majority  of  the  PAV  weight  budget.  These 
terms  are  inherently  coupled  and  more  investigation  is  needed  to  determine  PAV 
performance  limits  including  actuator  and  mechanism. 


Appendix  A 

Piezoelectric  Materials  and 
Actuators 

A.l  Piezoelectric  Actuators 

Piezoelectric  actuators  have  found  many  applications  in  macroscale  applications 
such  as  walking  microrobots  [81],  cooling  fans  [82],  inkjet  printers  [83],  and  mi¬ 
croscale  applications  such  as  RF  MEMS  switches  [84],  micropumps  [85],  and  MEMS 
gyroscopes  [86],  for  example.  Historically,  piezoelectric  materials  have  been  used 
for  ultrasonic  transducers,  resonators,  and  accelerometers.  This  section  gives  a 
brief  background  of  piezoelectric  materials,  and  describes  the  static  and  dynamic 
function  of  multilayer  beam  bending  actuators  and  novel  T-beam  actuators. 

A.  1.1  Piezoelectric  Materials 

Piezoelectric  materials  have  the  ability  to  convert  mechanical  energy  to  electri¬ 
cal  energy  (direct  effect)  and  vice  versa  (indirect)  and  are  available  in  many 
forms.  Quartz,  lead  zirconate  titanate  (PZT),  {l—x)Pb(Mgi/3Nb2/3)03—xPbTi03 
(PMN-PT)  and  polyvinylidene  diflouride  (PVDF)  are  examples  of  a  piezoelectric 
crystal,  ceramic,  relaxor-ferroelectric  and  polymer.  Figure  of  merits  for  these  piezo- 
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electric  materials  include  piezoelectric  coefficients,  electromechanical  coupling  co¬ 
efficient,  and  quality  factor.  The  linear  piezoelectric  constitutive  equations  relate 
the  electromechanical  properties  of  the  material 

Di  =  dijCTj  +  eJ^Ei  (A.l) 

Sj  =  sfjCTj  +  dijEi,  (A. 2) 


where  Sj  is  the  mechanical  strain,  <7*  is  the  mechanical  stress,  Ei  is  the  electric  held, 
Di  is  the  electrical  displacement,  sE  is  the  elastic  compliance  at  zero  electric  held, 
and  eJi  is  the  permittivity  at  zero  stress.  The  piezoelectric  coefficients,  are  third 
rank  tensors  which  in  reduced  tensor  notation  correspond  to  a  3  x  6  matrix[50].  In 
this  reduced  notation,  the  indices  (i  —  1  —  3)  define  normal  electric  displacement 
orientation,  ( j  =  1  —  3)  define  normal  mechanical  strains  and  (j  —  4  —  6)  represent 
shear  strains. 

The  electromechanical  coupling  coefficient  can  be  thought  of  as  energy  conver¬ 
sion  effectiveness,  and  is  defined  as  the  square  root  of  the  ratio  of  stored  mechani¬ 
cal  (electrical)  energy  over  input  electrical  (mechanical)  energy [83].  It  can  also  be 
written  as [52] 


k2  - 
Kij  ~ 


d2 

aij 

eTsE- 


(A.3) 


Mechanical  quality  factor  is  inversely  proportional  to  the  intensive  mechanical  loss, 
(tan S'),  of  the  material.  It  determines  the  sharpness  of  the  frequency  response 
of  a  piezoelectric  material  and  is  important  in  resonance  applications.  At  the 
macroscale,  bulk  crystals  and  ceramics  are  typically  used.  Microscale  sensors  and 
actuators  can  be  micromachined  from  bulk  material  or  deposited  using  thin  films 
(just  PZT). 

PZT  and  PMN-PT  are  considered  ferroelectric  materials,  which  have  a  sponta¬ 
neous  polarization  which  can  be  reversed  by  an  electric  field.  These  materials  have 
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a  Curie  temperature,  above  which  it  will  lose  its  polarization.  This  however  can 
be  reverse  by  re-poling  the  material  with  an  applied  held  while  near  or  above  its 
Curie  temperature,  and  allowing  the  material  to  cool.  The  material  can  also  lose  its 
polarization  due  to  an  applied  held  higher  than  its  coercive  held,  or  stress  which 
would  induce  this  held.  Ferroelectric  materials  are  useful  piezoelectrics  because 
the  polarization  can  be  designed  for  different  actuator  configurations,  however  the 
operating  temperature  must  be,  in  general,  at  least  one  half  of  the  Curie  temper¬ 
ature. 

PZT  has  seen  more  applications  because  of  its  ease  of  fabrication  either  in 
bulk  or  thin  him  form.  Despite  their  great  piezoelectric  properties,  PMN-PT  and 
other  relaxor  ferroelectric  materials  suffer  from  high  cost,  low  Young’s  modulus 
and  fracture  toughness.  PZT-4  and  PZT-5H  are  different  formulas  of  PZT  which 
can  be  considered  a  ’’hard”  or  ’’soft”  piezoelectric  material,  respectively.  In  general 
’’hard”  piezoelectrics  will  have  higher  quality  factor,  show  more  linearity,  and  be 
best  suited  for  resonance  applications.  ’’Soft”  piezoelectrics  on  the  other  hand, 
will  have  much  better  off-resonance  performance  but  show  significant  hysteresis  [50] 
and  much  lower  Curie  temperature.  PZT-5H  has  been  popular  in  NAV  actuator 
applications[44,  53]. 

The  indirect  piezoelectric  effect  is  used  for  piezoelectric  actuation.  If  there  is 
no  applied  stress,  Eq.  ??  can  be  simplified  to  A. 4 

€j  =  dijEi  (A. 4) 

where  6j  is  the  mechanical  strain  and  A,;  is  the  applied  electric  held.  Typically, 
the  index  3  is  used  to  indicate  the  direction  of  poling.  So,  an  applied  held  in  the 
direction  of  poling,  E3  will  cause  an  induced  strain,  e3  or  el5  proportional  to  the 
d33  and  d3 1  coefficients,  respectively.  Stack  actuators  use  the  d33  which  is  about 
three  times  larger  than  the  d31  and  is  capable  of  produce  large  forces  but  very 
small  displacement.  To  increase  displacement,  hextensional  and  cantilever  beam 
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actuators  are  used.  Flextensional  piezoelectric  actuators,  such  as  the  moonic, 
cymbal,  and  recent  uniflex[87]  use  an  elastic  element  to  amplify  displacement. 
These  actuators  use  the  d31  coefficient  to  amplify  tip  displacement  along  the  length 
of  the  beam.  The  following  section  gives  a  background  of  cantilever  beam  bending 
actuators. 

A.  1.2  Cantilever  Beam  Bending  Actuators 

Typical  PZT  cantilever  beam  bending  actuators  are  in  a  unimorph  or  bimorph 
configuration.  In  a  unimorph,  one  piezoelectrically  active  layer  bonds  to  another 
passive  layer.  An  applied  held  in  the  polarization  direction  will  cause  contraction 
in  the  active  layer  while  the  passive  layer  resists  the  deformation  causing  upward 
bending  in  the  beam.  Bimorphs  either  replace  the  passive  layer  with  a  piezoelec¬ 
trically  active  layer,  or  add  a  piezoelectric  layer  as  the  bottom  layer,  and  are 
configured  such  that  when  the  upper  layer  contracts,  the  bottom  layer  expands, 
increasing  the  performance  of  the  actuator.  Since  the  strain  induced  in  the  mate¬ 
rial  is  electric  held  induced,  the  voltage  required  to  actuate  is  dependent  on  the 
thickness  of  the  piezoelectric  layer.  To  reduce  this  operating  voltage,  multiple 
thin  layers  can  be  bonded  in  either  unimorph  or  bimorph  configuration  to  form 
multilayer  beam  bending  actuators[88]  cite  Ballas’s  book  too. 

Unimorph.  bimorph,  and  multilayer  beam  bending  actuators  have  been  studied 
extensively.  Models  of  the  static  and  dynamic  operation  will  be  presented  here  for 
completeness.  In  addition  to  the  geometry  of  the  piezoelectric  layer,  the  material 
choice  and  thickness  of  the  passive  layer  is  quite  important,  and  many  studies 
have  been  conducted  on  this  topic.  The  following  sections  describe  the  static  and 
dynamic  operation  of  bimorph  actuators. 

Free  tip  displacement,  blocking  force  at  the  tip,  and  stiffness  are  important  for 
the  static  operation  of  piezoelectric  actuators.  The  static  free  displacement  and 
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(c)  End  view 
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Figure  A. 1.1.  Schematic  of  T-beam  cantilever  beams:  (a)  top  view  showing  in-plane 
displacement  v(x,t)  and  tip  force  F%,  (b)  side  view  showing  out-of-plane  displacement 
w(x,t )  and  tip  force  ++  and  (c)  end  view  showing  cross  section  of  T-beam  with  design 
parameters  b,  web  width,  t.  flange  thickness,  s,  total  width,  and  h,  total  thickness 


blocking  force  equations  of  a  bimorph  with  center  shim  are  [74], 

3d31L2V  (5  +  1)  (25  +  1) 

2  (tpzt  +  tghimtpzt )  7l53  +  3 52  +  35  +  1 
3d3iEpztw(2tpzt  +  tshim)2V  25  +  1) 

8  Ltpzt  (5  +  l)2 


6f  = 

Fb  = 


(A.5) 

(A.6) 


where  A  =  Epzt/Eshim  and  5  =  tpzt/tshim ■  The  stiffness  of  a  a  bimorph  found 
as  k  —  Fs/Sf.  The  center  shim  is  required  to  be  electrically  conductive,  and  the 
bimorph  can  be  driven  in  multiple  configurations. 
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A.  1.3  T-beam  Actuators 

When  fabricating  cantilevers  using  PZT  film  deposition,  unwanted  curvature  due 
to  residual  film  stress  and/or  thermal  expansion  coefficient  mismatch  can  occur  [14, 
18].  For  certain  applications,  such  as  RF  switches  and  optical  alignment,  obtaining 
flat  cantilevers  is  crucial.  It  is  possible  to  compensate  for  the  curvature  during  the 
design  process,  however  it  requires  precise  knowledge  of  material  parameters  [19] 
which  can  be  difficult  in  practice.  A  more  robust  method  of  creating  flat  cantilevers 
is  by  increasing  the  substrate  thickness  to  withstand  the  residual  stress  at  the 
expense  of  tip  displacement [20].  This  however  leads  to  suboptimal  designs  for  the 
sake  of  a  flat  device.  Recent  advances  in  high  aspect  ratio  deep  reactive  ion  etching 
processes  [27-29]  have  made  possible  monolithic  fabrication  of  cantilever  T-beam 
actuators  from  bulk  PZT.  These  actuators  have  T-shaped  cross  section  that  can 
produce  out-of-plane  displacement  without  the  need  to  bond  or  deposit  multiple 
layers  [22], 

A. 1.3.1  Operation 

subsection  A.  1.3  shows  a  T-beam  cross  section  with  design  arameters  of  web  width, 
b ,  flange  thickness,  t,  total  width,  s,  and  total  height,  h.  Electrodes  are  deposited 
on  the  top  and  bottom  of  the  web  and  flange  regions.  The  T-beam  can  provide  in¬ 
plane  and  out-of-plane  displacement  by  selective  activation  of  the  web  and  flange 
electrode.  The  PZT  is  poled  through  the  thickness  from  top  to  bottom,  with  the 
bottom  electrode  acting  as  ground.  Assuming  uniform  electric  field,  a  positive 
electric  field  applied  to  the  web  in  the  direction  of  poling  causes  the  web  region  to 
expand  through  the  thickness  and  contract  longitudinally  by  the  dj,i  piezoelectric 
coefficient.  Similar  to  the  passive  layer  in  a  unimorph,  the  flange  regions  remain 
essentially  piezoclectrically  inactive,  constraining  the  lower  part  of  the  T-beam, 
forcing  the  beam  to  bend  upwards.  Similarly,  activating  both  flange  electrodes 
contracts  the  flange  region  while  the  web  region  resists  contraction  forcing  the  beam 
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to  deflect  downwards.  Activation  of  a  single  flange  causes  an  in-plane  deflection. 
Analytical  models  of  the  in-plane  and  out-of-plane  static  T-beam  tip  displacement 
in  response  to  applied  voltage  on  the  web  and  flange  electrodes  have  been  developed 
[89]  and  included  here  for  completeness.  The  in-plane,  Vf,  and  out-of-plane,  Wf, 
free  tip  displacement  in  response  to  an  applied  voltage,  V,  is  given  by 


vf 


wf 


Ap2d31L2V 

U2t 

_zEApd31L2V 

Eh,eff 


(A.7) 

(A.8) 


The  in-plane,  F2b ,  and  out-of-plane,  F3b ,  blocking  force  is  given  by 


F2  b  — 

F3b  = 


_ 3EAp2d3iV 

2 Lt 

2>zEApd3iV 

2 Z 


(A. 9) 
(A. 10) 


where  the  distance  from  the  top  of  the  flange  to  the  neutral  axis  is 

_  t2 s  -  t2b  +  2 tbh  -  bh2 
2  (ts  —  tb  +  bh) 

the  active  piezoelectric  cross  sectional  area,  Ap,  first  moment  of  area  about  the 
x2  axis,  Apo ,  area  moment  of  inertia  about  the  x2  axis,  I2l  and  effective  area  mo¬ 
ment  of  inertia  about  the  x3  axis,  /e//  are  given  in  table  1.  The  Young’s  modulus, 
E,  piezoelectric  strain  coefficient,  d3±,  permittivity,  633,  and  electromechanical  cou¬ 
pling  coefficient  of  PZT-4,  calculated  by  k3i  =  (Ed^/e^  ,  are  given  in  Table  A. 1.2. 
T-beam  cross  sections  have  been  optimized  for  displacement,  blocking  force  and 
mechanical  energy  output  [22],  Define  b*  =  b/s,  and  t*  =  t/s  as  the  web  and  flange 
ratio,  respectively.  Maximum  displacement  is  achieved  as  b*  and  t*  approach  0. 
There  exists  an  optimal  ratio  for  maximum  blocking  force,  b*  =  0.381,  t*  =  0.381, 
and  maximum  mechanical  energy  output,  b*  =  0.25,  t*  =  0.33.  This  will  be  useful 
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Table  A. 1.1.  Geometric  Parameters  of  T-beam 


Web 

Dual  flange 

Single  flange 

2e+h—  2t 

2  e—t 

2  e—t 

2h 

2 1 

2 1 

Ap 

bh 

t(s  —  b) 

\ t(s  —  b ) 

Ap2 

0 

0 

yis2  -  h2) 

ha 

A(s3-b3) 

b3h 

12 

A((2/i  -  t)b3  +  ts 3) 

Ip2 

b3h 

12 

Us3-V) 

^(s3-b3) 

Ib3 

( s—b)t 3 

12 

bh 3 

12 

+Ab(e-|)2 

+Afe(e  —  t  +  ())2 

-\-bh(e  —  t  T  ^ 

r 

bh 3 

(. s-b)t 3 

T  4-  bh3 

1b3 

12 

12 

1TJ3  +  12 

+A(e  —  t  +  |)2 

+Ab{e-lf 

— bh(c  —  t  -\~  y 

J3  ,eff 

(1  +  k^Ips  +  lb  3 

(1  +  k^lpz  +  Ib3 

(1  +  kl^IpZ  +  Ib3 

—kl1h2z2Ap 

—  k‘llt2z2Ap 

—  k21h2z2Ap 

in  designing  T-beam  actuators  for  applications  in  Nano  Air  Vehicles. 

A.  1.4  Static  and  Dynamic  Modeling  of  T-Beam  Actuators 

This  work  presents  static  and  dynamic  modeling  and  testing  of  PZT-5H  T-beam 
actuators.  Macroscale  T-beams  have  been  used  in  flapping  wing  mechanisms  driven 
statically [90]  and  at  resonance  [91].  The  results  of  this  work  can  also  be  applied 
to  microscale  T-beams  where  the  fabrication  process  is  particularly  advantageous. 
Dynamic  modeling  of  PZT-5H  T-beam  actuators  allows  the  T-beam  to  be  inte¬ 
grated  into  flapping  wing  mechanism  modeling  for  Nano  Air  Vehicle  (NAV)  appli¬ 
cations. 

A. 1.4.1  Fabrication  and  Experimental  Setup 

T-beam  actuators  are  fabricated  from  1  mm  thick,  finely  polished,  Cr/Au  elec- 
troded,  bulk  PZT-5H  using  a  high  precision  dicing  saw.  The  cross  section  of  the 
T-beam  has  web  width  b,  total  width  s,  flange  thickness  t,  and  total  thickness  h 
as  shown  in  subsection  A.  1.3(c).  The  material  properties  and  T-beam  dimensions 
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Table  A. 1.2.  PZT-5H  Material  Propertiesfl] 


Property 

PZT-4 

PZT-5H 

Young’s  Modulus,  E,  ( GPa ) 

78 

64 

Piezoelectric  Strain  Coefficient,  d3 1;  ( pm/V ) 

-122 

-265 

Relative  Permittivity,  633/60 

1700 

3400 

Density,  (kg/m3) 

7600 

7500 

are  shown  in  Table  A.  1.2.  After  the  T  cross  sections  are  realized,  the  individual 
T-beams  are  mounted  onto  a  glass  base  using  two  part  epoxy,  and  wires  are  sol¬ 
dered  to  the  web  and  bottom  electrode.  For  ease  of  handling  the  base  is  then 
bonded  to  4  mm  thick  FR4  PCB  board  material  which  can  be  mounted  onto  a 
Newport  linear  stage  as  shown  in  Figure  A.  1.2.  The  T-beams  are  poled  through 
the  thickness  of  the  web  region  using  2  V /  fxm  in  80  °C  peanut  oil  for  30  minutes 
and  then  cleaned  using  acetone,  isopropyl  alcohol,  and  detergent.  A  Model  609E-6 
high  voltage  Trek  amplifies  input  signals  from  a  Lab  VIEW  system.  A  Polytec  laser 
vibrometer  with  a  25  /im  diameter  laser  point  measures  static  tip  displacement  on 
the  T-beam.  An  Aurora  force  transducer  with  a  1  mm  diameter  glass  tube  mea¬ 
sures  blocking  force.  The  T-beam  is  clamped  to  a  Newport  linear  stage  using  the 
FR4  board  and  precisely  positioned  for  measurement  while  Lab  VIEW  records  all 
data.  Figure  ??  shows  a  photograph  of  the  mounted  beam,  the  laser  vibrometer 
used  for  displacement  measurement,  and  the  load  cell  used  for  force  measurement. 

A. 1.4. 2  Static  Results  and  Modeling 

During  the  blocking  force  measurement,  slight  flexibility  of  the  FR4  board  was  ob¬ 
served  using  the  laser  vibrometer.  This  softening  of  a  perfectly  clamped  boundary 
condition  can  be  modeled  as  a  free-free  beam  with  a  rotational  and  translational 
spring  at  the  base,  as  shown  in  Figure  A.  1.4.  The  Cartesian  coordinates  x  and 
z  denote  the  position  along  the  length  and  through  the  thickness  of  the  beam, 
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Wires 


Figure  A. 1.2.  Schematic  of  T-beam  mounted  onto  glass  base  and  FR4  board 


Table  A. 1.3.  T-beam  Model  Parameters 


Constant 

Expression 

a  2 

H1+ 

03 

(Ed3i)2(2e+h-2t)2 

8e33 

04 

Edsi  (2e-\-h—2t) 

2  h 

Os 

£33 

2  h2 

Tf , 

t(s  —  b) 

Ap 

bh 

h 

h 

yy  +  Ap  (e  -  t  +  §) 

Kv 

2^EApd31 

Jeff 

(h  +  «  +  ( I ,  -  Ap^^-) 

respectively.  Only  the  transverse  displacement  of  the  beam  will  be  considered  and 
is  denoted  by  w(x:t).  The  centroidal  axis  location  can  be  calculated  from  the 
T-beam  cross  section  parameters,  at  a  distance 

1  t2s  —  t2b  +  2  tbh  —  bh 2 


ts  —  tb  +  bh 


(A.ll) 
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Figure  A.  1.3.  Experimental  testing  setup  of  T-beam  actuator 
below  the  top  of  the  flange.  The  potential  energy  is 

U  =  Ub  +  Up  +  —kTw( 0,  t)~  +  —k rw  (0,  t )2  (A. 12) 

Ub=  f  f(- zw")2dV  (A.  13) 

Jvb  z 

Up=  f  HdV 
Jvp 


(A.  14) 
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Figure  A. 1.4.  Schematic  diagram  of  T-beam  with  modified  boundary  condition 

where  ()'  =  ^,  Ub  and  Up  are  the  potential  energy  of  the  passive  and  piezoelec- 
trically  active  region,  kx  is  the  translational  spring  constant,  is  the  rotational 
spring  constant,  and  E  is  the  Young’s  modulus  of  PZT-5H.  The  electric  enthalpy, 
H ,  is  given  as 


H  =  ( a2z 2  +  a3){w")2  +  aAVw"  +  a5V2.  (A.  15) 

where  the  coefficients  ai(i  =  2,  3, 4,  5)  are  given  in  Table  A. 1.3.  The  virtual  work 
done  by  the  force  generated  at  the  tip  is 


5W  =  FSw(L,t). 


(A. 16) 
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Substitution  of  the  potential  energy  into  the  principle  of  virtual  work,  f*f(6U  + 
SW)dt  =  0,  produces  the  field  equations 

w""  =  0,  Vx  G  (0,  L),  (A.  17) 


boundary  conditions 


EIeffw"(  0,  t)  +  KyV  —  kRw'(  0,  t)  —  0  (A.  18) 

EIeffw'"(0,t )  +  kTw'(0,t )  =  0  (A. 19) 

EIeffw"(L ,  f)  +  AyE  =  0  (A. 20) 

EIeffw"\L,t)  -  F  =  0  (A. 21) 


where  the  effective  bending  stiffness  of  the  T-bearn  EIeff  and  the  coefficient  Ky 
is  given  in  Table  A. 1.3.  The  solution  for  Equation  A. 17  is 

w{x)  =  c3x3  +  c2x2  +  C\ x  +  c0.  (A. 22) 


The  coefficients  q(z  =  0, 1,2,3)  are  solved  using  the  boundary  conditions  in  Eqs. 
(??)  -  (??).  The  free  displacement  and  blocking  force  at  the  tip  of  the  beam  are 


wf 


Fb 


KyL2V 
2  EIeff 

3KvE2V 

2  A3  +  6  EIeff(L2/kR  +  1  /hr) 


(A. 23) 
(A. 24) 


The  d3 1  coefficient  given  by  the  manufacturer  in  Table  A.  1.2  is  measured  at  low 
electric  fields.  Nonlinear  electromechanical  effects  begin  to  contribute  significantly 
when  soft  PZT  materials  are  actuated  at  fields  greater  than  0.005  V j fim[ 92].  The 
piezoelectric  strain  coefficient,  d33  can  be  considered  a  function  of  electric  held, 
temperature,  and  stress[93]  and  can  increase  40%  in  response  to  high  electric 
fields [94].  In  this  work,  the  focus  is  to  develop  useful  models  for  T-bearn  actu- 
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ators  over  a  certain  operating  regime.  Figure  A.  1.5  shows  the  experimental  and 
theoretical  peak  values  of  static  displacement  using  a  40%  larger  d3i  than  the  data 
from  the  manufacturer  in  Table  A. 1.2.  Figure  A. 1.6  shows  the  experimental  and 
theoretical  peak  value  of  the  blocking  force,  where  the  parameters  k r  and  k t  were 
adjusted  to  simulate  the  deflection  of  the  FR4  board  thereby  reducing  the  force 
predicted  in  the  model. 


Figure  A. 1.5.  Experimental  (circle)  and  theoretical  (solid)  tip  displacement  versus 
electric  field 
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Figure  A.  1.6.  Experimental  (circle)  and  theoretical  (solid)  blocking  force  versus 
electric  field 

A. 1.4. 3  Dynamic  Modeling 

The  kinetic  energy  of  the  T-beam,  neglecting  rotational  inertia,  is 

1  CL 

T  —  -  pAw2dx.  (A. 25) 

2  Jo 

where  p  is  the  density  of  PZT-5H,  and  A  =  Ab  +  Ap  is  the  total  cross  sectional 
area  of  the  beam.  Using  the  potential  energy  term  derived  in  (Equation  A.  14), 
the  virtual  work  (Equation  A.  16),  and  the  kinetic  energy  (Equation  A. 25)  into 


104 


Hamilton’s  principle  f*f(SW  +  5T  —  5U)dt  yields  the  equations  of  motions 

EIeffw""  +  pAw  =  0  (A. 26) 

with  the  same  boundary  conditions,  Eqs.  (A.  19)  -  (A. 21).  The  equations  of  motion 
can  be  modified  to  include  strain-rate  damping  and  viscous  air  damping  by  [95] 

EIeffw""  +  csw""  +  caw  +  pAw  =  0  (A.  27) 

where  cs  is  the  strain-rate  damping  coefficient,  ca  is  the  viscous  air  damping  co¬ 
efficient.  Considering  the  small  area  of  the  T-bearn,  assume  the  air  damping  will 
have  neglegible  effect  compared  with  the  strain-rate  damping  (ca  =  0).  Taking  the 
Laplace  transform  of  equation  (??)  with  ca  =  0  yields 

pAs2w  +  cssw""  +  EIeffw""  =  0.  (A. 28) 


Assume  a  solution 

w(x,  s )  =  C\  cos (/3x)  +  C2  sin (f3x)  +  C3  cosh(^x)  +  C4  sinh(/3x)  (A. 29) 


where 


P(s) 


(  s2pA  y/4 

V  css  +  EIeff) 


(A. 30) 


and  s  =  juj.  Solving  for  the  coefficents  Ci(i  =  0, 1,2,3)  using  the  boundary 
conditions  (A.19-A.21)  yields  the  solution  for  tip  displacment  in  the  frequency 
domain 


N(P) 

D(P) 


V(s) 


w(L,  s ) 


(A. 31) 
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where 


N(P)  =  Kv  {cos  {PL)  -  cosh  (PL))  (A.32) 

\  kRkR 

+  - — (sinli  (PL)  —  sin  (PL))  (A. 33) 

kri 

+(w+^b)sin(/JL)8inh(/Ji)))  (A-34) 

+  (  - - h  J  cosh  (PL)  sin  (fiL)  (A. 35) 

\kn  kR  J 

+  (  - ]  cos  (PL)  sinh(/3L) )  (A. 36) 

\kR  kR  J  J 

D(P)  =  P(  1  +  cos  (PL)  cosh  (PL))  (A. 37) 

+  ^  (cos  (PL)  cosh  (PL)  —  1)  (A. 38) 

(A. 39) 

\  kR  kT  J 

x  (cosh(/3L)  sin  (PL)  +  cos  (PL)  sinh  (PL))  (A. 40) 


Figure  A. 1.7  shows  the  experimental  and  theoretical  frequency  response  of  the 
T-bearn  actuator.  The  damping  coefficient  cs  was  adjusted  to  match  the  exper¬ 
imental  data.  During  strong  AC  electric  fields,  soft  piezoelectrics  also  exhibit 
softening  due  to  elastic  nonlinearity [92],  The  Young’s  modulus  was  reduced  by 
10%  from  Table  A.  1.2  to  simulate  this  effect. 

A. 1.4. 4  Conclusion 

This  work  presents  static  and  dynamic  testing  and  modeling  of  PZT-5H  T-bearn 
actuators.  These  models  are  useful  tools  in  further  studies  involving  NAV  flapping 
wing  mechanisms.  The  results  include  boundary  condition  modifications  to  simu¬ 
late  imperfect  clamping.  This  could  be  useful  because  the  mechanical  grounding 
in  a  lightweight  NAV  airframe  may  not  provide  a  perfect  clamped  boundary  con¬ 
dition.  Nonlinear  effects  on  the  electromechanical  properties  of  soft  PZT  under 
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Figure  A. 1.7.  Experimental  (solids)  and  theoretical  (dashed)  displacement  versus 
frequency  in  response  to  0.01  V / /jm  electric  field 

strong  fields  were  considered,  and  piezoelectric  coefficients  and  Young’s  modulus 
were  modified  based  on  experimental  data.  A  frequency  response  solution  has 
been  derived  including  strain-rate  damping.  The  accuracy  of  the  models  greatly 
depend  on  proper  experimental  data  to  fit  key  parameters,  but  is  useful  for  future 
modeling  efforts. 


Appendix  B 

Design  and  Extra  Results  of  the 
Lionfly 

B.l  Single  Degree  of  Freedom  Design 

B.1.1  Kinematics 

The  following  section  studies  the  effect  of  the  mechanism  parameters  on  static 
operation  of  the  LionFly  mechanism.  The  device  designed  to  flap  at  a  maximum  of 
90  deg.  at  resonance,  which  sets  the  flexure  lengths,  La,  Lb,  Lo  for  a  given  flexure 
thickness,  tf.  The  starting  position  of  the  actuator,  xref,  is  chosen  such  that  the 
starting  flapping  angle  is  45  deg.  The  flexure  lengths  and  mechanism  parameters: 
link  lengths,  Lab,  Loa,  and  vertical  offset,  r4,  must  satisfy  Equation  4.1.  In  the 
first  step  in  the  design  optimization  it  is  reasonable  to  assume  a  certain  actuator 
displacement,  because  flapping  angle  is  based  purely  on  kinematics.  Later,  we 
include  the  spring  forces  from  the  flexures.  As  the  width  of  the  actuator  increases, 
the  actuator  bending  stiffness  and  blocking  force  production  increases  and  the  two 
results  should  converge.  First,  understanding  the  effect  of  Lab,  Loa,  and  vertical 
offset,  r4,  on  flapping  angle  is  made.  For  Figure  B.1.1  and  Figure  B.1.2,  assume 
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P-P  Flapping  Angle  Varying  r  with  ±  100  Lt  m  x  ^  (deg.) 


Figure  B.1.1.  Peak  to  peak  flapping  angle  versus  crank  length.  Top  plot,  Lab  = 
1.81  mm,  bottom  plot  r 4  =  1  mm. 


+/-  100  fim  displacement  from  the  actuator  (0.05%  strain  of  a  15  mm  bimorph  or 
applied  150  V). 

Figure  B.1.1  shows  the  peak  to  peak  flapping  angle  as  a  function  of  crank  length 
and  varying  vertical  offset  and  coupler  link  length.  In  the  first  subplot  using  r 4, 
a  coupler  link  value  was  chosen  that  leads  to  least  amplification  Lab  —  1.81  mm, 
and  the  same  idea  for  the  coupler  link  plot,  holding  74  =  1  mm.  The  coupler 
length  doesn’t  greatly  affect  the  flapping  angle  magnitude.  Setting  the  crank 
length  to  zero  (r2  is  nonzero  because  of  the  flexure  lengths),  will  yield  maximum 
amplification. 

In  Figure  B.1.2  we  see  that  although  the  peak-peak  flapping  angle  increases 
when  v 4  goes  to  zero,  it  becomes  more  nonlinear.  Increasing  r 4  shifts  the  dead 
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Figure  B.1.2.  Flapping  response  with  prescribed  bimorph  displacement,  centered  on 
xref  such  that  <fi  =  45  deg..  Notice  how  the  =  0  response  begins  to  curve 


position  of  the  mechanism,  when  the  crank  length  and  coupler  length  become 
colinear,  to  a  negative  angle.  So,  a  nearly  linear  response  is  achieved  when  we 
make  74  large  (1  mm  is  easy  to  fabricate  using  a  thick  glass  slide).  The  responses 
are  centered  around  45  deg.,  and  have  different  xref  positions.  The  coupler  to 
crank  length  ratio,  Lab / L<ja  determines  the  deviation  from  a  perfect  sinusoid  the 
output  will  have.  The  higher  the  ratio,  the  closer  to  a  pure  sinusoid  the  output 
will  be. 

B.1.2  Static  Response 

Now  we  solve  the  equilibrium  equation  for  the  actual  bimorph  position,  using 
Equation  4.9.  The  coupler  length  and  vertical  offset  were  fixed  at  2  mm,  and 
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Figure  B.1.3.  4>pp  versus  crank  length,  Loa  with  applied  voltage  of  ±150  V,  cj) o  =  45°: 
(solid)  calculated  i f>  by  solved  bimorph  displacement,  (dashed)  4>  calculated  by  bimorph 
free  displacement 


1  mm  respectively.  The  constraints  for  the  length  of  the  crank  link,  Loa ,  actuator, 
Lact,  and  width  of  the  actuator,  wact  are 


0  <  Loa  <  1  mm 

(B.l) 

5  <  Lact  <  25  mm 

(B.2) 

0.5  <  wact  <  10  mm. 

(B.3) 

The  flexure  thickness  is  tf  =  25/xm,  widths  Wo  =  6  mm,  wa  =  wb  =  3  mm,  and 
Ef  =  2.8  GPa.  With  and  applied  voltage  ±150  V,  Equation  4.9  is  solved  for  x, 
and  Equation  4.2  is  used  to  find  0(x). 

The  results  show  the  effect  of  the  flexure  stiffness  in  reducing  peak  to  peak  flap- 
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Width  =  10.00  mm 


Figure  B.1.4.  cj)pp  versus  actuator  length,  Loa  with  applied  voltage  of  ±150  V,  </>o  = 
45°:  (solid)  calculated  <fi  by  solved  bimorph  displacement,  (dashed)  <j)  calculated  by 
bimorph  free  displacement 


ping  angles  for  most  actuator  geometry  and  crank  lengths.  The  results  show  that 
for  any  given  crank  length,  the  wider  actuators  will  approach  the  free  displacement 
values  as  expected.  For  narrow  actuators,  smaller  crank  lengths  (increased  kine¬ 
matic  amplification)  do  not  buy  larger  flapping  angles.  For  a  given  width,  there  is 
a  critical  length  after  which  the  flapping  angle  decreases,  as  shown  in  Figure  B.1.4. 
With  maximum  width,  maximum  length,  and  shortest  crank  length  the  largest 
flapping  angles  are  achieved,  however,  even  then  the  free  displacement  values  are 
not  reached  (for  tf  —  25  /im). 
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Figure  B.1.5.  <j)pp  versus  actuator  width,  Loa  with  applied  voltage  of  ±150  V,  (j) o  = 
45°:  (solid)  calculated  <fi  by  solved  bimorph  displacement,  (dashed)  (j)  calculated  by 
bimorph  free  displacement 

B.1.3  Amplitude  Dependent  Linear  Aerodynamic  Drag  Force 

The  aerodynamic  drag  force  on  a  single  degree  of  freedom  flapping  wing  mech¬ 
anism  is  a  quadratic  nonlinear  damping  term  dependent  on  the  flapping  veloc¬ 
ity.  Developing  a  linear  approximation  of  this  term  will  allow  design  of  flapping 
wing  mechanisms  based  on  resonance  operation  and  include  amplitude  dependent 
damping.  The  device  has  maximum  flapping  amplitude  and  angular  velocity  at 
resonance,  however  it  is  difficult  to  estimate  what  the  magnitude  is  without  ade¬ 
quate  experimental  results  to  tune  damping  coefficients.  Following  the  approach 
of  Davis  [96],  we  can  approximate  a  linear  damping  coefficient  using  the  formula 
for  the  aerodynamic  drag  force  on  the  wing,  as  a  function  of  the  forcing,  linear 
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Figure  B.1.6.  Contour  plot  of  (j>DC  amplitude  over  actuator  length  and  width. 


damping  and  natural  frequency.  If  we  focus  on  the  single  degree  of  freedom  system 
with  flapping,  0,  and  temporarily  neglect  -0,  the  kinetic  energy  is  simply 

T  =  — ( Jyy  +  Js  +  Jam)<p2  —  (B.4) 

where,  Jyy  is  the  inertia  about  the  flapping  axis,  Js  is  the  inertia  of  the  wing  spar 
about  the  flapping  axis,  and 

Jam  =  I  jc{r)2r'2dr.  (B.5) 

The  virtual  work  is 


5VF  =  GactVSx  +  FNSr  CPx  +  c^Scj), 


(B.6) 
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where  is  the  linear  viscous  damping  in  vacuum,  and  using  Lagrange’s  equations 
with  <f>  as  the  generalized  coordinate,  we  can  write  the  equations  of  motion  as 


where 


0  +  2  Hl4>  +  2/Xq|0|0  +  UJn  —  Fq(V) 


PL  =  wf  = 

ZJ(j) 


rcpxpairCN(g )  f0R  c(r)r'2dr 
2( Jyy  +  Js) 


PQ  ~ 


K  =  -2-  I  kact  +  -d-fs(x) 


dx‘ 


X=XQ  , 


a  — 


d(j){x) 


dx 


X=XQ 


(B.7) 


(B.8) 

(B.9) 

(B.10) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 


From  the  results  in  [96],  if  the  forcing  function  F0  is  sinusoidal  at  the  resonant 
frequency  con,  we  can  arrive  at  an  estimate  for  the  amplitude  at  resonance 


<b(V) 


-3flL7T  +  y(3  jlL^)2  +  48/1qFo(F)7T 
16nQvjn 


(B.15) 


from  which  we  can  calculate  an  voltage  dependent  linear  damping  coefficient, 


(q(V) 


8[iQar(V) 
37 r 


(B.16) 
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such  that  the  equations  of  motion  can  be  rewritten  as 

4>  +  2  (Cl  +  C Q{V))j>  +  co2n  =  F0(V)  (B.17) 

To  validate  this  approach,  LF07  was  tested  using  125  V DC  and  multiple  voltage 
amplitudes  (25  VAC,  50  VAC,  and  75  VAC)  in  air  and  in  vacuum.  ??  shows  the 
flapping  frequency  response  taken  using  high  definition  photographic  measurement 
and  image  processing.  The  results  show  consistent  prediction  of  resonant  amplitude 
as  voltage  levels  and  flapping  amplitudes  are  increased,  validating  the  approach 
to  estimate  amplitude  at  resonance.  This  can  now  be  used  to  estimate  flapping 
velocity  at  resonance  as  a  function  of  wing  shape,  resonant  frequency,  and  forcing. 

B.1.4  Rotational  Hinge  Stiffness 

To  achieve  lift,  we  must  have  rotation  in  proper  phase  with  flapping.  This  requires 
properly  designed  flexure  hinges  that  allow  passive  rotation  in  response  to  aerody¬ 
namic  forces  acting  upon  the  wing.  To  predict  this  force,  the  steady  state  angular 
velocity  of  the  wing  must  be  found.  Using  the  results  of  last  section,  assuming  a 
sinusoidal  response,  the  angular  velocity  can  be  written  as 

<t>  =  unap.  (B.18) 

If  the  wing  is  extremely  light,  the  force  balance  of  the  system  in  Fig.  ??  is 

rcpyFN  -  k =  0  (B.19) 

Substituting  Eq.  ??  and  ??  into  ??,  we  can  rearrange  terms  to  solve  for  the  hinge 
stiffness  for  a  desired  ip 

u  -  rcpypN 


(B.20) 


Figure  B.1.7.  Photographed,  image  processed  and  theoretical  flapping,  <p(t)  versus  fre¬ 
quency  at  (125  VDC/  25  VAC)  (left),  (125  VDC/  75  VAC)  (middle),  and  (125  VDC/ 
75  VAC)  (right):  Experimental  4>(t)  in  air  (blue  circles)  and  vacuum  (green  diamonds), 
and  theoretical  <f>(t)  in  air  (blue  solid)  and  vacuum  (green  solid). 


The  length  of  the  rotational  flexure  hinge  is  determined  by  the  maximum  bending 
angle  equation  ref  TBD.  Since  the  thickness  is  limited  by  fabrication  constraints, 
we  control  the  stiffness  of  the  hinge  using  the  width.  To  achieve  the  desired  k^,, 
the  flexure  width  must  be 


^  flex 


12k^Lfiex 

7  KeEfie:ct3flex 


(B.21) 
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B.2  Additional  LionFly  Data 

New  LionFly  designs  were  fabricated  to  minimize  flexure  width  according  to  Equa¬ 
tion  B.21,  and  with  an  extra  photolithography  step  to  define  a  intermediate  flexure 
thickness  level  for  vein  definition.  Yield  of  devices  during  fabrication  drastically 
decreases  when  attempting  to  create  flexure  thickness  below  20  /j,m.  Addition¬ 
ally,  successfully  fabricated  devices  were  extremely  fragile  and  had  very  short  life 
cycles,  preventing  sufficient  testing.  Even  at  these  reduced  thicknesses,  testing 
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Figure  B.2.1.  Combined  photograph  of  Tipulidae  1  at  rest  and  actuated  at  -150  VDC. 
Thin  flexure  hinges,  a  75  deg.  initial  angle,  and  25  mm  long  bimorph  contributed  to  an 
increase  in  static  performce  with  33  deg.  flapping  for  -150  V. 


showed  nearly  in  phase  flapping  and  rotation,  showing  that  the  ’’‘light  wing”’  as¬ 
sumption  is  invalid.  Figures  B.2.1  and  B.2.2  show  the  static  and  dynamic  response 
of  the  LionFly  modeled  after  a  Cranefly,  Tipulidae ,  an  insect  of  the  order  Diptera. 
Figure  B.2.3  shows  the  static  response  of  the  LionFly  modeled  after  a  Hoverfly, 
Eristalis,  also  of  the  order  Diptera.  This  device  shows  the  greatest  static  perfor¬ 
mance  of  any  LionFly  to  date.  Finally,  a  scaled  down  version  of  LF1203  with  10 
mm  span  was  fabricated  and  tested  for  lift.  The  device  shows  a  force  transducer 
measured  110  fiN,  which  corrected  for  leverage  is  approximately  60  /j,N,  showing 
slightly  less  lift  than  the  LF  1203. 
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Figure  B.2.2.  Photograph  of  flapping  at  75  V  and  45  Hz  for  Tipulidae  1.  Large 
amplitude  flapping  at  75  deg.  demonstrates  the  effectiveness  in  flapping  mechanism  per¬ 
formance.  Improvements  can  be  made  to  extend  the  cycle  life  of  these  devices,  however. 


B.3  Review  of  SU-8  Material  Properties 

SUEX  is  a  relatively  new  material  and  little  literature  exists  that  report  material 
property  testing  as  a  function  of  processing  conditions.  SU-8  negative  photoresist 
shares  many  similar  characteristics  to  SUEX  and  has  numerous  articles  that  report 
the  influence  of  processing  conditions  on  material  properties.  The  key  material 
properties  we  are  interested  in  are  Young’s  modulus  and  yield  stress.  Table  B.3.1 
shows  a  list  of  SU-8  material  properties,  and  associated  fabrication  process  notes. 
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Figure  B.2.3.  Combined  photograph  of  Eristalis  1  at  rest  and  actuated  at  -150  to  150 
VDC.  Thin  flexure  hinges,  a  80  deg.  initial  angle,  and  25  mm  long  bimorph  contributed 
to  an  increase  in  static  performce  with  71  deg.  total  flapping  angle,  the  best  recorded 
for  any  LionFly. 
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Lift  Frequency  Response  in  Air,  100  V,  100  Vdc-  LF1218R10 


Figure  B.2.4.  Average  value  lift  versus  frequency  for  LF1218R10  with  10  mm  wingspan 
actuated  at  (100  VDC /  100 VAC)  over  three  independent  runs. 
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